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EXECUTIVE  SUMMARY 


OBJECTIVES 

Describe  hem  the  field  strength  of  an  HP  signal,  expressed  in  decibels  (dB)  above  or  below 
1  microvolt  per  mettr  reference,  is  calculated  in  different  HP  propagation  prediction  programs 
and  how  the  accuracy  of  the  predicted  field  strength  values  frenn  these  programs  can  be  deter¬ 
mined  and  presented. 

RESULTS 

All  of  these  prediction  programs  produce  median  predictions  of  the  rms  field  strength.  Seven 
HP  propagation  prediction  programs  have  been  reviewed  for  the  Polar,  Equatmial,  Near  vertical 
incidence  E]q)eriment  (PENEX)  project.  These  include  three  empirical  based  programs  (Medusa 
PROPHET,  FTZ,  and  PTZ4),  a^  four  analytical  programs  (HFTDA,  lONCAP.  ASAPS,  and 
AMBCOM).  AMBCOM  is  the  cmly  ray  tracing  program  induded. 

The  implementation  of  a  data  screening  program  DASCR3  allows  the  devdopment  and 
gmieration  of  a  powerful  statistical  description  of  die  characteristics  of  die  measured  field 
strength  and  of  how  well  the  seven  candi^te  inograms  predid  observations.  It  offers  all  the 
statistical  requirements  suggested  by  the  Comite  Consultadf  International  des  Radiocommunica¬ 
tions  (CCIR)  for  the  determinatiiMi  of  the  accuracy  of  a  field  strength  prediction  program.  (The 
CQR  is  the  International  Radio  Consultation  Committee  in  English.)  Useful  statistical  parame¬ 
ters  produced  by  DASCR3  that  can  describe  the  accuracy  of  a  predicted  field  strength  value 
indude  die  following:  average  residual  (bias),  root-mean-square  residual  (or  standard  deviation), 
average  relative  residual  (relative  bias),  root-mean-square  relative  residual,  average  absolute 
relative  residual  (magnitude  of  the  error  in  die  model),  ccMrelation  coefficient  between  observed 
and  predicted  values,  standard  enxH'  of  the  estimate  of  linear  regression,  and  the  constants 
necessary  to  rqiresent  the  residual  distribution  by  a  Johnson  distribution  and  its  conespmuling 
test  of  fit  informatiem.  DASCR3’s  ability  to  allow  and  st(»re  up  to  40  different  auxiliary  variables 
allows  the  comparison  to  be  subdivided  into  many  subcategories.  DASCR3  usage  also  allows  the 
determination  ^  possible  improvements  diat  might  be  made  to  these  field  strength  prediction 
IMTograms. 

RECOMMENDATIONS 

1.  The  accuracy  of  the  field  strength  predictions  from  the  seven  HP  predictions  described 
hmein  be  determined  using  PENEX  data. 

2.  For  comparison  to  other  accuracy  detotninations,  these  same  seven  programs  be 
compared  to  CCIR  Data  Base  D.l  using  DASCR3. 

3.  The  accuracy  of  these  programs  be  detomined  as  a  function  of  the  following  auxiliary 
variables:  the  month;  year;  sunspot  number;  circuit  identifier  (name  and  path  transmitter  and 
receiver  coordinates);  frequency;  great-circle  distance;  24  predicted  hourly  values  (always 
monthly  median  values)  of  sky-wave  field  strength  in  dB  relative  to  1  pV/m,  of  path  basic  MUF, 
of  tile  percentage  of  the  days  po*  month  when  the  frequency  is  below  Ae  path  basic  MUF,  of 
solar  zenith  angle  and  cosine  of  the  solar  zenith  angle  at  path  mi(4>oint,  of  E-layer  MUF;  24  pre¬ 
dicted  hourly  values  at  each  reflection  point  (control  point)  of  E-layo'  critical  frequency,  secant 
of  tiie  angle  of  incidence  on  the  D-layer,  critical  frequency  of  the  sporadic-E  layer  foEs,  and  the 
sporadic-E  layer  blanketing  frequency  fbEs;  and  the  eight  3-hour  magnetic  index  Kp  values. 


4.  The  Jdmscm  distributicm  parameters  be  determined  for  each  value  of  auxiliary  variable 
used  to  determiim  the  accuracy  of  each  inx)9tun. 

5.  The  results  of  die  accuracy  stwlies  of  these  programs  be  used  to  make  recmnmenda* 
tima  oa  hofw  each  program  might  be  im[H'Oved. 
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INTRODUCTION 


This  report  describes  how  the  field  strength  of  an  HF  signal,  expressed  in  dB  above  or  below 
a  1  microvolt  per  meter  reference,  is  derived  in  different  HF  prediction  pro^ams,  particularly 
the  ones  that  will  be  evaluated  in  Project  PENEX.  HF  sky-wave  field  strength  is  a  measure  of 
the  signal  intensity  as  it  appears  at  the  receiver  antenna.  Field  strength  is  one  of  the  least 
uiuierstood  parameters  characterizing  HF  signals  and  is  one  of  the  most  difficult  parameters  to 
measure  by  calibration  techniques.  The  report  will  present  a  general  discussion  of  HF  prediction 
jnrograms  and  their  differences;  it  will  also  provide  a  detailed  discussion  of  how  each  program 
develops  their  field  strength  numbers.  A  substantial  amount  of  text  is  devoted  to  the  develop¬ 
ment  of  the  basic  transmission  loss  equations  used  in  each  program  to  show  the  difference 
between  the  models. 

The  final  section  will  jivesent  how  the  observed  and  predicted  data  will  be  compared  in  the 
analytical  studies  and  will  describe  a  special  data  screening  program  used  for  that  purpose  at  the 
Naval  Command,  Control  and  Ocean  Surveillance  Center,  Research,  Development,  Test  and 
Evaluation  Division  (NRaD). 

HF  PREDICTION  TECHNIQUES 

The  following  section  is  taken  from  AGARDograph  No.  326  (pp.  69-72, 1990).  To  deter¬ 
mine  the  performance  of  an  ionospheric-dependent  radio  system,  more  than  just  an  ionospheric 
model  is  needed.  The  ionospheric  model  must  be  tied  to  a  set  of  equations  or  a  formulation  that 
enables  the  simulation  of  the  propagation  of  radio  waves  through  the  ionospheric  model.  The  set 
of  equations,  or  the  formulation  chosen,  together  with  the  ionosj^eric  model,  are  often  termed 
an  ionospheric  propagation  model.  When  the  ionospheric  model  that  is  contained  in  the 
propagation  m^el  can  be  used  for  making  predictions  of  the  ionospheric  structure,  the  propaga¬ 
tion  model  is  termed  an  ionospheric  propagation  prediction  model.  The  propagation  model  must 
provide  the  method  for  calculating  the  geometry  pertinent  to  the  radio  system,  as  well  as 
methods  for  handling  information  about  required  performance  levels:  transmitter  power,  signal 
level  and  modulation,  antenna  characteristics,  receiver  location,  and  noise  environment.  The 
degree  that  each  of  these  is  incorporated  into  the  propagation  model  often  determines  the 
complexity  of  the  model. 

Most  of  the  HF  propagation  models  available  assume  that  signals  are  reflected  from  the 
ionospheric  E-  and  F-regions  according  to  strict  geometrical  considerations.  The  ionospheric 
parameters  at  the  reflection  points  are  estimated  from  the  ionospheric  model  and  are  u^  as 
input  to  the  formulation  relating  to  the  reflection  process.  The  details  of  the  method  used  to 
evaluate  the  reflection  of  HF  signals  from  the  ionospheric  regions  (i.e.,  the  evaluation  of  modes) 
vary  with  different  propagation  models. 

Performance  predictions  are  made  for  many  purposes,  such  as  system  design,  frequency 
management,  and  operational  improvements.  Most  of  the  propagation  methods  were  originally 
intended  to  provide  information  of  a  long-term  predictive  nature  by  using  monthly  median 
predictions  of  ionospheric  structure;  however,  a  trend  has  emerged  in  recent  years  to  utilize 
propagation  predictions  on  much  shorter  time  scales.  The  complexity  of  the  long-term  and 
short-term  jvopagation  prediction  methods  is  generally  as  different  as  the  approaches  used. 

The  best  known  long-term  performance  prediction  methods  involve  the  used  of  large- 
scale  computer  programs.  The  work  of  Lucas  and  Haydon  (1966)  was  the  first  long-term. 
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computer-based  i^ogram  of  its  smt.  The  concepts  of  stfvice  probability  and  reliability  were 
introduced  in  this  fu’ogram,  HFMUFS,  aiul  was  subsequently  rq>laced  by  that  of  Bai]^usai  et 
al.  (1969).  Thoe  were  four  distinct  versitms  of  these  programs,  dubbed  ITS78.  Each  was  given  a 
slightly  different  name  and  was  color  coded  (Red  De^  Blue  Deck,  Yellow  Deck,  Buff  Deck) 
according  to  the  cokr  of  the  cards  on  which  it  was  sent  out. 

The  final  version,  HFMUFES4,  (Haydon  et  al.,  1976)  gained  international  usage.  Since  ite 
introduction,  the  lONCAP  program  (Teters  et  al.,  1983)  has  also  become  widely  used.  These 
jM'Ograms  provide  the  means  to  calculate  HP  iMrq>agation  parameters  m  any  location  on  the  earth. 
Held  strength,  mode  reliability,  and  Maximum  Usable  Fr^uency  (MUF)  are  but  a  few  of  the 
parametors  that  are  obtained  from  these  programs.  They  enable  the  program  user  to  specify 
antenna  gains  as  a  function  of  take-off  angle  and  to  specify  required  systems  perfntnance,  in 
terms  of  the  signal-to-noise  ratio  evaluated  at  the  receiving  point  of  the  circuit.  Both  programs 
have  common  features,  such  as  use  of  the  same  sets  of  numerical  coefficients  to  represent  the 
morphological  behavitST  of  the  itmospheric  structure  and  the  atmospheric  noise  eiq^ected  at  the 
reception  point.  There  are,  however,  significant  differences  among  HFMUFS,  HFMUFES4,  and 
lONCAP.  The  majtu*  changes  from  HFMUFS  in  HFMUFES4  are  as  follows: 

1 .  All  numerical  coefficients  rqiresenting  the  ionospheric  characteristics  were  calculated 
as  functions  of  universal  time. 

2.  E-layer  propagation  characteristics  woe  calculated  frmn  numerical  coefficients  repre¬ 
senting  E-layer  critical  frequencies  (Leftin,  1976). 

3.  Numerical  coefficients,  rqiresenting  the  minimum  virtual  height  of  the  F-region, 
were  included  for  calculating  the  semi-thickness  of  the  F-layer  (Lefdn,  Ostrow  & 
Preston,  1967). 

4.  Revised  values  of  manmade  noise  and  its  frequency  dqiendence  were  included. 

5.  A  method  fot  combining  two,  or  more,  noise  sources  of  nearly  equal  amplitudes  was 
added. 

6.  A  new  fcnmula  fix'  estimating  absoptitm,  including  a  winter  anomaly  effect,  was 
derived  (Schultz  &  Gallet,  1970). 

7.  The  chi-square  distribution  was  used  to  evaluate  all  distributions  Q^acharisen  & 

Crow,  1970). 

8.  Revised  excess  system  losses  were  inclined. 

9.  System  perfcxmance  predictions  could  be  made  for  sporadic  E-prq>agation. 

10.  The  numerical  miq)s  of  foP2  were  continuous  in  month  and  sunspot  number. 

11.  Numerical  coefficients  rq)resenting  atmospheric  noise,  as  a  function  of  universal 
time,  woe  included  (Zacharisen  &  Jones,  1970). 

12.  Numerical  maps  of  the  continents,  for  use  in  ground  loss  calculations,  wo-e  added 
(Zacharisen,  1972). 

13.  Provision  was  made  to  use  up  to  three  different  transmitting  and  receiving  anteimas 
over  the  HF  band. 

14.  Modifications  were  made  to  allow  antenna  patterns  to  be  read  into  the  program. 

The  iont^pheric  loss  term  in  lONCAP  differs  significantly  from  that  used  in  either  HFMUFS 
or  HFMUFES4.  The  lONCAP  uses  the  same  set  of  ft,F2  coefficients  as  does  the  HFMUFS. 
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The  original  vision  lONCAP  used  the  same  1-MHz  representation  of  atmospheric  noise.  The 
omrent  lONCAP  uses  a  1-MHz  rqjresentation  of  atinosph^c  noise  due  to  Spaulding  and 
Washburn  (1985).  The  varicnis  mo(tes  in  these  three  programs  are  computed  in  different 
manners;  ccmsequently,  the  signal-to-noise  ratio  that  is  calculated  by  each  of  the  programs  for 
the  circuit  conditions  is  different.  The  lONCAP  program  has  a  distinct  advantage  of  tire 
Barghausen  et  al.  (1969)  prp^am  by  enabling  the  user  to  incorporate  into  the  calculation 
q)ecific  knowledge  about  the  ionosphoe,  such  as,  an  electron  density  profile  obtained  from 
independent  infmmation. 

In  utilizing  a  prc^agation  prediction  method,  the  user  must  specify  the  particulars  of  the 
circuit,  such  as  the  transmitter  and  receiver  location,  transmitter  power,  transmitter  and  receiver 
anteima,  and  the  quality  df  so^ce  required.  In  addition,  the  universal  time,  month,  and  sunspot 
numbor  that  are  apprq^riate  for  the  period  i<x  which  calculations  are  to  be  perfcumed  must 
specified.  There  are  numerous  ouq>ut  q)tions  that  are  available,  including  Maximiun  Usable 
Frequency  (MUF)  for  the  circuit,  l^e  Lowest  Useful  Frequency  (LUF),  and  the  field  strength  fm* 
any  fiequency  that  has  been  indicated  by  the  user.  The  mode,  signal-to-noise  ratio,  predicted 
signal  reliability,  and  take-off  angle  for  each  mode  are  likewise  available. 

The  three  prediction  programs  discussed  above  are  complete  HP  propagation  performance 
prediction  programs.  There  is  an  existing  class  of  programs  that  can  considered  a  subset  of 
these.  These  programs  are  concerned  primarily  with  evaluating  the  field  strength  of  an  iono¬ 
spheric-dependent  radio  system.  Models  of  this  type  are  given  in  CCIR  repcxt  252-2  (CCIR, 
1970a),  the  supplement  to  CCIR  Rqxxt  252-2  (CCIR,  1980)  and  CCIR  Rqxxt  894  (CCIR, 
1982);  CCIR  !^p<nt  894-2  (CCIR,  1990a).  The  field  strength  calculations  given  in  Rq>ort  252-2 
(CCK,  1970a)  are  consistent  with  the  method  used  in  HFMUFS.  On  the  other  hand,  the  field 
strength  calculations  that  are  given  in  the  Supplement  to  Report  252-2  (CCTR,  1980)  are  more 
complex  than  the  method  of  calculation  of  field  strength  used  in  ITS78  or  HFMUFES4.  The 
complexity  is  due  to  a  significant  difference  in  the  manner  in  which  the  ionospheric  modes  are 
evaluated.  In  the  CCIR  Report  252-2  approach  to  ionospheric  reflection  estimations  and  mode 
evaluation,  the  pertinent  calculations  are  perfixmed  at  specific  points,  called  control  points, 
along  the  propagation  path  determined  by  the  path  lengA.  No  account  is  taken  of  the  change,  (x 
gradient,  in  electron  d^ity  along,  or  transverse  to,  the  propagation  path.  These  gradients  are 
accounted  fix  in  the  Supplement  to  Reptxt  252-2.  It  thus  provides  a  more  physically  ^pealing 
calculation  at  an  increase  in  computational  time  by  a  facUx  of  10  to  30. 

The  field  strength  predictirm  method  given  in  CQR  Repcxt  894  had  as  its  roots:  wrxk 
poformed  by  CCIR  Interim  Wtxking  Party  6/12  to  develop  a  sky-wave  propagation  prediction 
program  fix  use  in  planning  the  HF  broadcasting  service  (TTU,  1984).  This  field  strength  model 
is  actually  a  combination  of  two  field  strength  programs;  a  simplified  version  of  CCIR  Rqxxt 
252-2  is  used  for  path  lengths  of  less  that  7(XX}  km,  and  the  field  strength  model  developed  by 
Deutsche  Bundespost  (FTZ)  (Damboldt,  1976)  is  used  for  distances  of  greater  than  9000  km.  A 
linear  interpolation  scheme  is  employed  fix  distances  between  70(X)  and  9000  km.  M(xe  recent 
revisions  to  this  method  are  desoib^  in  CCIR  Reptxts  894-1  (CCIR,  1986)  and  894-2  (CCIR, 
1990). 

F(x  the  propagation  models  given  above,  the  field  strength  is  evaluated  fix  each  mode  that  is 
determined  according  to  the  geometry  incorpcxated  into  the  program.  The  selection  of  the  modes 
that  are  chosen  to  determine  the  overall  field  strength  f(x  a  given  frequency  is  not  the  same  for 
each  of  the  programs.  Generally,  however,  three  or  four  of  the  modes  that  are  associated  with  the 
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least  amount  of  loss  are  chosen;  then,  the  antenna  gain  is  incorporated  into  the  field  strength 
calculation  for  each  mode. 

The  FTZ  propagation  program  (Damboldt,  1976;  Damboldt  &  Suessmann,  1989)  employs  an 
empirical  field  strength  calculation  that  is  based  upon  observations  collected  over  a  number  of 
HF  circuits,  most  of  which  terminate  in  Germany.  The  data  that  have  been  gathered  for  more 
than  10  years  for  certain  circuits  are  related  to  predicted  ionospheric  critical  frequencies  to 
obtain  an  empirical  method  for  determining  Held  strength.  In  particular,  the  field  strength 
recordings  revealed  a  steady  increase  in  signal  strength  from  the  LUF  to  a  maximum  value, 
following  approximately  an  inverse  frequency  dependence.  This  frequency  is  called  the  LUF  in 
this  method;  a  formula  for  its  calculation  is  provided.  It  is  different  than  the  classical  LUF 
calculation.  After  the  maximum  value  is  reached,  the  field  strength  decreases  until  it  reaches  the 
operational  MUF,  which  is  higher  than  the  classical  MUF.  This  is  the  consequence  of  several 
mechanisms  that  are  not  taken  into  account  by  the  other  prediction  techniques.  Because  of  its 
simplicity,  the  FTZ  method  was  adopted  in  the  mid-1970s  for  propagation  prediction  programs 
designed  for  computers  with  limited  memory  in  the  early  desktop  models,  such  as  PROPHET 
and  early  versions  of  Medusa.  The  weakness  in  FTZ  models  is  that  the  peak  of  the  field  is 
dependent  on  how  the  operational  MUF  and  LUF  are  used  to  determine  the  field  strength.  If 
these  frequencies  are  well  chosen,  then  the  field  strength  is  for  the  minimum  hop  mode,  and  the 
antenna  gain  can  be  determined  for  it.  Otherwise,  the  field  strength  {^diction  will  be  inaccurate 
at  any  given  frequency,  and  the  mode,  for  which  the  field  strength  is  represented,  will  be 
unknown.  How  the  MUF  and  LUF  are  chosen  for  the  FTZ  model  will  be  discussed  in  some 
detail  later. 

Report  894  (CQR,  i982)  provides  the  basis  of  yet  another  {vopagation  model  that  was 
developed  at  the  First  Session  of  the  HF  Broadcasting  Conference  (ITU,  1984).  This  model, 
referred  to  as  the  HFBC84,  was  developed  specifically  for  planning  the  use  of  the  HF  spectrum 
for  broadcasting  proposes.  The  primary  difference  between  HFBC84  and  the  Report  894  model 
is  in  the  manner  that  the  antemia  gain  is  taken  into  account  in  the  computation  of  field  strength. 
Before  the  selection  of  the  modes,  which  are  to  be  combined  to  determine  field  strength  of  a 
given  frequency  on  paths  of  less  than  7000  km,  the  anteima  gain  at  the  appropriate  take-off  angle 
for  each  mode  is  added  to  the  field  strength.  The  resultant  field  strength  is  determined  by  using 
the  strongest  E-mode  and  the  two  strongest  F-modes  for  paths  up  to  4000  km.  Between  4000  and 
7000  km,  only  the  two  strongest  F-modes  are  considered.  For  paths  greater  than  9000  km,  the 
maximum  anteima  gain  that  occurs  between  0-  and  8-degrees  elevation  angle  is  used  in  the  field 
strength  computation.  The  inclusion  of  the  anteima  gain  in  the  field  strength  calculation,  {nior  to 
the  selection  of  the  modes  that  are  combined  to  form  the  resultant  field  strength,  leads  to  a  much 
improved  field  strength  prediction.  The  HFBC84  program  provides  an  efficient  means  to 
determine  the  area  serviced  by  the  HF  broadcast  transmitter  and  to  assess  the  likely  interference. 

Yet  another  prediction  program  was  produced  by  the  International  Working  Part  (IWP)  6/1 
of  the  CCIR  in  response  to  Recommendation  No.  514  (HFBC-87),  which  invites  the  CQR  to 
“...undertake  studies  of  the  propagation  prediction  method  adopted  by  the  Conference  and  to 
recommend  both  improvements  in  the  method  and  later,  if  necessary,  an  improved  method  to  be 
used  in  the  future.”  An  interim  report  was  produced  by  the  CQR  as  Report  894-2  (CQR,  1990). 
In  1991,  a  final  report  (CQR,  1991)  presented  eight  elements  of  the  HraC-87  propagation 
method  for  which  improvements  were  recommended. 

There  is  a  class  of  prediction  programs  that  differs  considerably  from  those  discussed  in 
previous  paragraphs,  lliese  programs  are  concerned  primarily  with  tracing  the  signal  rays 
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through  the  ionosphere.  The  one .  gram  of  this  type  to  be  considered  in  the  PENEX  evalua¬ 
tions  is  called  Ambient  Ionospheric  Communication  Predictions  at  HP  Program  AMBCOM 
(Hatfield,  1980;  Hatfield  et  al.,  1987;  Smith  &  Hatfield,  1987).  AMBCOM  was  bom  out  of  work 
done  by  Stanford  Research  Institute  in  the  1960s  to  develq)  computer  codes,  called  NUCOM 
Codes,  that  would  simulate  ionospheric  changes  after  above  ground  nuclear  explosions. 

In  AMBCOM,  the  ionosphere  is  modeled  in  three  parabolic  layers.  Ionospheric  tilts  and 
initial  frequency  gradients  are  taken  into  account  by  specifying  the  parabolic  parameters  at  as 
many  as  41  points  along  the  path.  These  parameters  were  initially  derived  from  the  Institute  for 
Telecommunications  Sciences  coefficients  used  in  HFMUFES.  Then  they  were  modified  to 
incorporate  a  high-latitude  ionospheric  model  (Elkins  &  Rush,  1973a  and  1973b;  Vondrak  et  al., 
1978),  an  auroral  absorption  model  (Vondrak  et  al.,  1978),  and  a  sporadic-E  model  ( Phillqjs, 
1963;  Sinno,  Kam,  &  Kirukawa,  1976;  Kolawole,  1978).  If  desired,  actual  measurements  may 
be  used  in  place  of  parameters.  The  propagation  analysis  consists  of  a  rt^id,  semi-analytic, 
two-dimensional,  ray  tracing  routine  based  on  the  Kifi-Fooks  method  (Westover  &  Roben, 

1963).  Both  topside  and  bottomside  reflections  from  the  normal  ionospheric  layers  are  allowed. 
AMBCOM  computes  propagation  losses  with  a  homing  feature  for  evaluation  of  specific 
point-to-point  communication  circuits,  along  with  binary  error  rates  and  signal-to-noise  ratio. 

EMPIRICAL  HF  FIELD  STRENGTH  MODELING— 1976  TO  PRESENT 

INTRODUCTION 

In  1976,  the  original  development  of  the  MINIMUF  model  (Levine,  1976)  provided  simple 
formulation  to  calculate  the  HF  MUF  using  minicomputer  and  emerging  microcomputer 
technology.  Empirical  HF  oblique  sounder  data  provided  the  basis  of  the  MINIMUF  and  QLOF, 
a  simple  LUF  model.  At  the  same  time,  work  in  Germany  by  Deutsche  Bundespost  produced  a 
simplified  empirical  field  strength  (Damboldt,  1976).  This  model,  called  Nachdchtentechische 
Zeitschrift  (FIZ),  used  data  that,  for  the  most  part,  had  been  collected  over  long  paths  (i.e., 
greater  than  7000  km).  Because  the  model  was  simple  and  easy  to  use,  this  last  fact  was  just 
ignored.  Given  that  the  LUF  and  MUF  boundaries  could  be  estimated,  it  was  generally  felt  that 
FTZ  gave  a  reasonable  estimate  of  predicted  field  strength.  The  field  strength  model,  coupled 
with  a  noise  model  that  was  develqied  later,  led  to  predictions  of  signal-to-noise.  In  1976,  the 
first  PROPHET  system  was  devel(q)ed  on  an  ANAiYK-3  militarized  minicomputer  that  had 
32-Kbytes  (KB)  of  RAM.  Over  the  years,  literally  dozens  of  versions  of  PROPHET  were 
developed  by  using  the  MINIMUF,  QLOF  and  FTZ  models  as  the  heart  of  the  predictions. 
Traditional  methods  to  calculate  these  parameters  were  just  too  cumbersome;  however,  computer 
technology  over  the  last  several  years  has  made  that  a  moot  point.  The  question  of  how  good 
these  predictions  were  has  lingered  throughout  this  period,  and  the  issue  has  largely  remained 
unresolved. 

It  should  be  remembered  that  in  1976,  MINIMUF,  the  genesis  of  these  simple  models,  was 
developed,  such  that  MUF  calculations  could  be  done  in  less  that  32  KB  of  RAM.  The  original 
MINIMUF  was  80  Basic  statements  in  length.  The  ^proach  is  untraditional  because  MINIMUF 
is  an  emulation  on  how  the  MUF  boundary  fluctuates  as  a  function  of  time  of  day,  season,  and 
sunspot  cycle.  Simply  put,  it  is  an  empirically  calibrated,  dual  RC  lag  circuit.  Its  developers  were 
engineos.  As  a  result  of  scientists  in  the  HF  propagation  field  taking  every  opportunity  to 
discredit  MINIMUF,  the  accuracy  of  MINIMUF  itself  has  been  well  documented  by  its 
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developing  organization  (Sailors,  Moision,  &  Brown,  1981;  Sailor,  Sprague  &  Rix,  1986;  Roy 
&  Sailors,  1987).  The  same  reaction  was  not  true  for  QLOF  because  the  LUF  boundary  is  a 
fuzzy,  vague  function  that  is  signal-to-noise  dependent.  No  one  could  come  up  with  a  better  idea. 
Last,  FTZ  was  devel(^)ed  to  stq)port  the  long  range  shortwave  te-oadcasting  service  of  Deutsche 
Bund^post.  The  data  that  it  used  in  its  empirical  develq)ment  all  had  one  common  feature:  one 
end  point  was  in  Germany.  This  immediately  opens  the  question  as  to  its  global  applicability. 
Even  with  its  detractrvs,  HIOPHET  continues  to  be  one  of  the  favored  tools  by  certain  military 
users,  some  commercial  broadcasters,  and  the  amateur  shmtwave  radio  community. 

Subsequent  sections  will  discuss  the  evolution  of  FTZ;  its  use  in  PROPHET, which  will  be 
referred  to  as  FS;  the  expansion  to  HFTDA  as  a  successor  model  in  PROPHET  and  Medusa;  and 
the  new  improved  FTZ4  field  strength  models.  In  addition,  the  mme  traditional  programs, 
lONCAP,  AS  APS,  and  the  raytrace  program  AMBCOM  will  be  discussed. 

THE  FTZ  MODEL 

The  basic  field  strength  calculation  developed  in  the  original  version  of  FTZ  in  1976  remains 
the  same  today.  It  was  first  derived  by  Beckmaim  in  1965  (Beckmaim,  1965, 1967)  and 
described  the  variation  in  field  strength  within  a  transmission  range  as 


E  = 


-30  +  Gg  +  10  log  p 


(1) 


where 

^  =  Sky-wave  ims  field  strength  in  dB  above  1  pV/m 

=  free  space  field  strength 
=  Gain  of  the  transmitting  antenna 
P  =  Effective  radiated  power  (eip) 

=  Operational  MUF,  or  the  upper  frequency  limit  for  a  transmitter 
power  of  1  Mw  erp,  and  a  receiving  field  strength  of  1  pV/m  (or 
1  Kw  erp  and  30  dB  below  1  pV/m) 

=  frequency  where  the  field  strength  of  a  1-Mw  transmitter  is 
1  pV/m. 

FTZ’s  q>proach  to  field  strength  calculation  u  a  traditional  q)proach  by  using  CCIR  formulation 
with  some  “tweaking”  in  the  and-^m  calculations  by  using  locally  acquired  empirical  data.  The 
mistake  most  people  make  by  using  FTZ  and  FTZ4  is  to  assume  thati^j  is  the  classical  LUF,  and 
that-ffl  is  the  classical  MUF.  In  fact,  they  are  not  and  are  quite  different. 

Damboldt  (1976)  stressed  that  the  Beckmann  fcxmula  empirically  comprises  all  different 
factors  influencing  propagation.  It  yields  only  an  estimated  field  strength.  The  modes,  solar 
zenith  angle,  angle  of  incidence,  blanketing  and  other  phenomena,  are  not  dealt  with  separately, 
as  is  the  case  in  the  mm-e  analytical  field-strengdi  prediction  methods.  These  phenomena  are 
contained  partly  in  the  characterization  of  the  circuit  frequency  boundaries  with  fi  and  fg ,  the 
classical  MUF  used  in  the  determination  of  .  The  empirically  determined  K-factor  of  the 
MUF  computation  implicitly  comprises  all  other  influences  on  the  field-strength  calculation. 
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One  of  the  major  strengths  of  empirically  derived  models  is  that  the  data,  used  for  the  modeling, 
inherently  contain  the  subtle  variations  that  are  so  difficult  to  model  explicitly.  In  many  cases, 
empirical  models  will  out-perform  a  model  that  uses  more  traditional  analytical  approaches. 

Damboldt  (1976)  observed  that  as  the  operating  frequency  increased  away  from  the  fre¬ 
quency  where  the  signal  is  no  longer  useful,  the  field  strength  increased  to  a  point  under  the 
MOF,  and  then  started  decreasing  past  the  classical  MUF  to  a  frequency  called  the  operational 
MUF  (OMUF).  The  causes  of  this  controversial  higher  boundary  included  scatter,  transequato- 
rial  propagation,  field-aligned  irregularities,  sporadic-E,  off-great-circle  propagation,  and 
ducting;  all  phenomena  that  are  not  easily  modeled,  but  known  to  exist.  The  OMUF  is  derived 
by 


(2) 


with 


1  +  ff- 


fg,  noon 


+  X 


fg,  noon 


(3) 


where 

s  Classical  MUF 
s  gyrofrequency 

fg,  noon  s  fg  for  local  noon  of  the  respective  control  point. 
fg,  min  =  minimum  hourly  value  of  the  classical  MUF. 

The  constants  W,X,Y  are  chosen  empirically;  they  are  dependent  on  the  geographical  position  of  the 
HF  path.  To  a  certain  degree,  they  take  into  account  propagation  through  the  auroral  zones  and  the 
geomagnetic  equator.  Values  are  shown  in  table  I  in  Uie  original  Damboldt  paper  (Damboldt,  1976) 
and  are  given  here  in  table  1. 

Table  1.  Empirical  values  of  W,X,  and  Y  constants  used  for  the 
determination  of  the  correction  factor  K. 


Orientation 

Constants 

W 

X 

Y 

East-West 

0.1 

1.2 

0.6 

North-South 

0.2 

0.2 

0.4 

The  theory  behind  by  using  the  K-factor  is  briefly  presented.  The  dynamics  of  the  ionosphere 
are  characterized  by  diurnal  and  seasonal  variations  of  the  classical  MUF.  On  the  other  hand, 
ionospheric  irregularities  become  more  apparent  as  the  F-layer  height  increases  and  its  critical 
frequency  decreases;  consequently,  the  highest  value  for  K  is  obtained  in  winter  nights  when  the 
classical  MUF  values  are  at  their  lowest. 
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The  variability  of  the  first  two  terms  in  the  formula  of  the  K-facttx  is  due  to  the  ratio  of  the 
hourly  values  of  Ae  standard  MUF  to  the  noon  value  of  MUF  and  its  reciprocal  value.  re^)ec' 
lively.  The  term  with  the  coefficient  X  describes  the  increase  of  the  K-fzctx  from  day  to  night; 
however,  since  during  die  day,  the  X-facttx’  has  to  reach  a  certain  value,  the  term  with  the 
coefficient  W  is  added,  counteracting  the  term  described  frrst.  That  last  term,  with  the  coefficient 
Y,  varies  with  die  ratio  of  the  minimum  standard  MUFf^,  min  to  the  noon  value  fg ,  noon,  and 
it  expresses  the  increase  of  the  K-f actor  from  winter  to  summer. 

The  classical  MUF  is  based  on  the  CCIR  atlas  of  ionospheric  characteristics  (CCIR,  1970b). 
The  Control  Point  Method  is  qiplied.  In  the  case  of  the  E*  and  F-layers,  the  two  control  points 
are  respectively  1000-  and  2000-kxn  away  from  the  terminals.  After  the  great-circle  path  and  the 
great-drcle  distances  to  the  control  points  have  bemi  determined,  the  24-hourly  values  of  fg  are 
calculated  according  to  the  CCIR  atlas. 

The  calculation  of  fi  is  dmved  from  the  formula  for  non-deviadve  absorption,  which  is  the 
basis  of  the  CCIR  LUF  calculation.  The  parameter  fi  is  found  by  setting  the  free-space  field 
strength  for  1000  kW  eip,  nomalized  to  1  kW,  minus  the  non-deviadve  absorption  fen*  the 
minimum  hop  mode  loss  equal  to  30  dB  below  1  |xV/m  and  is  given  by 

52  Vcosxd  +  0.009*Ri2) 

/  =  5  3  j  .M. _ 


where 

I 

a  seasonal  factor  that  is  also  dependent  on  the  geogr^hical  posi¬ 
tions  of  both  terminals  and  can  be  determined  from  /  =  fj 
(National  Bureau  of  Standards,  1948) 

X 

= 

solar  zenitii  angle 

Ri2 

= 

the  12-month  running  mean  of  sun^t  number 

the  angle  of  incidence  at  the  D-layer 

Dp 

the  oblique  path  length 

the  number  of  penetration  points 

= 

the  gyrofrequency. 

With  the  above  formula  for  f  i  ,  the  LUF  is  detoinined  for  the  daylight  hours.  During  the  night,  the 
LUF  is  assumed  to  be  dq)endent  only  on  the  distance  between  transmitter  and  receiver  so  that 


An  = 


\  3000 


(5) 


In  addition,  as  there  is  a  certain  lag  between  the  time  of  sunset  and  the  decrease  of  D-region  ioniza¬ 
tion,  the  decay  from  day-LUF  to  night-LUF  is  accounted  for  the  three  hours  after  sunset  by 

fi  =  (6) 
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With  the  determination  of  andfj ,  the  field  strength  can  be  calculated  by  calculating  the 
free-space  fudd  strmgth  using 


=  20  log  3x10^^  (7) 

wheieP  istheerp.and  X>  is  the  distance  between  transmitter  and  receiver.  This  is  the  field  strength 
produced  by  1-kW  ii^ut  to  a  short  dipole  over  perfect  ground.  The  gain  of  the  transmit  antenna  is 
4.8  dB,  relative  to  an  isotrqjic.  To  determine  die  field  strength  relative  to  an  isotropic,  it  is  necessary 
to  subtract  4.8  dB  from  the  field  strength  given  by  equation  (1). 

THE  FTZ4  MODEL 

In  1989,  a  revision  to  FTZ,  called  PTZ4  (Damboldt  &  Suessmann,  1989)  was  introduced.  It 
follows  the  same  tqiproach  as  FTZ,  in  that  it  uses  the  formula  in  equatitm  (1),  requiring  the 
definition  of  fj  and  .  There  were  some  modifications  in  how  these  functitms  were  derived. 
The  parameter  K  in  equation  (3)  now  has  a  lead  ctmstant  of  1.2  instead  of  1.  The  value  of  f  j  for 
the  three  hours  after  sunset  is  now  given  by 

f,  =  (8) 

The  majo'  changes  in  FIZ4  include  an  iffl[X‘oved  MUF  derivadtm  model,  consideration  of 
E-Regitm  intervention,  and  some  correction  factors  to  consider  different  ranges. 

FTZ4  introduces  the  FTZMUP2  MUF  prediction  model.  It  derives  f  ( 3  0 0  0 ) ,  for  a 

given  location  from  updated  CQR  data  bases,  and  calculates  a  MUF(4000),  as  a  function  of 
time,  by  the  eiqiression 


MUF  -  1.1  ♦  f^2  *  «(3000)  (9) 

The  foP2  and  3  000 )  values  are  inteqiolated  as  a  function  of  season,  local-time,  geomagnetic 
latitude  for  fo^2 ,  and  geographic  latitude  for  ^  ( 3 0 0 0 )  fitun  tables  (teveloped  from  CCIR  Rq>ort 
430-4  Atlas  of  Ionospheric  Characteristics  (CGR,  1983).  The  accura^  of  ^  ^^roach  is  a  direct 
function  on  the  size  and  accuracy  of  the  tables  of  coefficients.  As  is  the  case  in  Ae  presentation  of 
any  of  these  mini-HF  prediction  {n'ograms,  there  is  the  comparison  to  MINIMUF  to  demonstrate  its 
supoicnity — and  this  paper  is  no  excq>ti(m.  It  is  not  clear  whether  or  not  they  used  the  MINIMUF 
(H*  MINIMUF  3  J  for  the  comparison;  however,  they  claim  an  improvement  by  reducing  the  standard 
deviation  in  FTZMUF2  by  a  factor  of  two.  It  is  a  moot  point,  as  MINIMUF  3.5  has  been  rq>laced 
by  MUF85,  which  has  demonstrated  an  even  larger  improvement. 

A  second  improvmnent  in  FIZ4  is  the  introduction  of  the  E-Region  MUF.  By  using  a 
function  determined  by  Rawer  (1952, 1956),  the  E(D)MUF  is  calculated.  The  critical  frequen¬ 
ces  (E)  for  the  E-layer  (in  MHz)  are  determined  as  follows: 

fg  =  Kg*coa^x 

Kg  =  2.25  +  l.Scos^  +  (0.01  -  0.007  COS^)  (10) 

n  =  0.21  +  0.12COS^  +  0.0002R^2 
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where 


^12  -  12-m(mth  runniqg  mean  of  the  sunqxjt  number 

X  s  sun’s  zenith  angle 

^  s  geographical  latitude. 

The  24*hourly  values  of  the  E-layer  E(D)MUF  are  determined  frtmi 

E(D)MUF  (11) 


with  the  distance  factm 

Eo  ^  {  ( (  (  { <-4.368*10-9*£rarc+1.335*10-‘^)  *  Ear  c-S  .911  *10-^) 

*£arc+0 . 0002625)  •Earc-0. 005039)  *Earc+0 . 03761)  (12) 

♦Earc-0 . 01332)  *Eazc-*^0 .2065 


^i^iere  Bare  =  E-layer  hop  length  in  radians. 

The  basic  MUF  for  the  whole  circuit  is  chosen  by  taking  the  higher  value  of  the  E-  or 
F2-MUF  for  each  control  point,  and  then  taking  the  lowest  value  of  all  three  control  points. 

The  last  change  in  FIZ4  is  the  addition  of  aconection  factor  to  the  qrerational  MUF  (f.) 
that  is  dqrendent  cm  path  lengtti.  It  has  been  known  for  a  long  time  that  the  Damboldt  database  is 
predmninantly  long  paths.  In  the  new  model; 


f  =  E  *  f  y  ♦  C, 


(13) 


where  C,  equals  2-(D/4000)^,  for  circuits  <4000 km,  or  1,  for  circuits  >4000  km.  This  has  the  effect 
of  reducing  the  operatiotud  MUF  for  short  paths. 

A  careful  check  of  all  the  other  constants  used  throughout  the  FTZ4  indicated  that  it  uses  the 
same  constants  as  wo-e  used  in  the  oiginal  FTZ,  excq>t  that  the  lead  constant  in  equatitm  (3)  is 
now  1.2  instead  of  1.0. 

FTZ4  comes  on  a  stand-alone  3.S-inch  floppy  and  is  suitable  for  use  on  a  standard  desktop 
mini-computer.  In  the  program  comparisons  for  Project  PENEX,  the  FTZ4  model  will  be  used. 
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THE  MEDUSA  PROPHET  HP  FIELD  STRENGTH  PREDICTION  MODELS 


The  calculation  of  field  strengdi  in  die  PROPHET  HP  Signal  Assessment  Systems  has  been 
under  continual  change  since  1976  as  i»w  and  imjx'oved  af^voadies  were  devel(q)ed. 

Tlw  oiginal  PROPHET  syston  was  developed  to  operate  in  a  cmnputer  system  that 
ccmtained  a  very  limited  amount  of  memoiy,  eidto*  32  KB  or  64  KB  RAM;  therefore,  when 
FTZ  first  af^ieared,  the  differential  between  die  classical  MUF  and  the  operatimi  MUF  was 
viewed  wiA  some  skepticism.  The  first  versions  of  PROPHET  took  a  sinqile  q^oach,  they 
used  the  MUF  calculated  widi  MINIMUF  as  the  and  the  LUF  value  calculated  by  die 
original  QLOF  model  as  f  j .  Over  time,  it  was  learned  diat  whm  ctnnpared  with  small  ammmts 
field  strength  data,  the  peak  in  the  calculated  FS  was  occurring  too  low  in  the  usable  band¬ 
width  and  the  adoptiem  of  an  operational  was  necessary. 

A  simplified  method  of  estimating  die  operational  MUF  was  implemented  (Sailors,  1990).  It 
was  assumed  that  the  MUF  value,  at  a  given  time,  seasem,  geogr^hic  position,  and  sunspot 
number,  is  a  statistical  parameter  and  has  a  Gaussian  distribution.  Given  this  assumptiem,  the 
operadmial  MUF  in  PROPHET  was  determined  by  the  jnroduct  of  the  MUF  and  the  99.1 
percentile  value  of  the  MUF  distributiem  as  exjxessed  by 


f.  =  1.85  *  HPF  •  MUF 


(14) 


where  HPF  is  the  90di-percentile  value  factor  for  an  assumed  Gaussian  MUF  distribution  (Baighau- 
sen  et  al.,  1969).  The  reason  that  this  equatkm  was  used,  instead  of  the  approach  used  in  FTZ,  is  that 
the  K-factor  was  based  mainly  <hi  data  measured  on  paths  terminating  at  one  site  (Germany),  and 
might  contain  aqiects  specific  to  that  site.  Table  2  lists  the  HPF  (Fg)  and  die  FOT  (Fj)  as  a  function 
of  tune  of  day. 

As  revised  versions  of  MINIMUF  were  introduced,  they  were  adopted  into  the  PROPHET 
FS  calculadmi.  Between  1976  and  the  present,  the  models  used  were  MINIMUF,  M1N1MUF-3.S 
and  MUP85  (Sailors  et  al.,  1986),  the  latter  bdng  the  model  presently  in  die  Medusa  {Hoject 
PROPHET.  MUF85  rqvesented  a  significant  inqirovement  over  M1NIMUF-3J5,  especially  at 
high  latitudes  (Sailors  et  al.,  1986). 

The  value  of  the  LUF,  f  j  ,  was  derived  ftxim  the  anpirical  QLOF  series  of  models  (Sailors  & 
Moision,  1987),  the  latest  of  these  being  QLOF  2.0,  which  is  used  in  Medusa  PROPHCT.  Excqit 
for  these  revisions,  the  basic  equation  (1),  as  shown  in  the  FTZ  discussion,  is  the  equation  that  is 
IM^reently  being  used  in  Medusa  PROPHET.  Fm-  fi«quencies  above  the  MUF,  the  fidd  strength  is 
reduced  by  a  loss  of  12  dB  po-  MHz  until  a  maidmum  increase  in  loss  of  60  dB  is  obtained. 

HFTDA~A  HYBRID  APPROACH  TO  FIELD  STRENGTH  PREDICmON 

The  mnpirical  approaches  of  the  1970s  and  the  mid-1980s  were  driven  by  the  need  for 
simplicity  to  accommodate  limited  computer  meminry.  Since  1987,  desktop  computers  have 
become  faster  and  have  mov  memoy,  so  the  miginal  requirement  fra-  simplidty  has  been 
overcome  by  events.  A  new  qiproach  was  adcqited  at  the  Naval  Command,  Control  and  Ocean 
Surveillance  Center  (NCCOSC),  Research,  Develqiment  and  Test  Division  (NRaD),  formerly 
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Tible  2.  Factors  for  calculating  POT  and  HPP  from  tbe  MUF  and  the  local  time. 
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die  Naval  Ocean  Systems  Center  (NOSC),  to  develt^  a  new  field  strengdi  model  to  rqilace  the  w- 
sims  of  FTZ,  which  had  been  in  use  for  about  15  years,  for  cme  that  would  more  Kcurately  charac- 
toize  amditions  at  path  lengths  of  less  than  7000  km  (Sailors,  1990;  Systems  Exploradmi,  Inc., 
1990).  The  new  mo^  called  i^FIBA,  revoted  to  more  classical  methods  by  using  system  loss 
equatUms  to  calculate  the  predicted  path  loss  between  tnnsadttar  and  receiver.  The  E^layer,  P-layer, 
and  mixed  mode  hops  are  calculated.  For  path  lengths  greato-  than  1,000  km,  the  FTZ  model  tued 
in  PROPHET/It^usa  was  retained.  At  ttese  l<mger  ranges,  QLOF  Vosion  2.0  is  used  for  obtaining 
fi ,  and  fm  is  found  by  using  MUFs  produced  by  MINIMUF-8S  in  equation  (14). 

HPTDA  ccmtains  a  model  relying  (m  {xopagation,  via  the  regular  E-layer,  as  was  first 
described  in  ESSA-ITSAl  (Lucas  &  Haydon,  1966)  and  the  F2  layo*  that  was  develqied  in  the 
semi-enq>irical  modeling  fn*  MINIMUF-8S  (Sailors  et  al.,1986).  The  program  finds  the  lowest 
order  mode  for  the  E-layer,  die  P-layer,  and  a  mixed  E-layer  and  P2-layer  mode.  Only  erne 
E-layer  hop  is  ermsidered  in  the  mixed  mode.  The  remaining  hops  are  via  the  P2-layer. 

The  E-layer  MUF  is  calculated  as  follows:  first,  the  ionospheric  absorption  index  is  deter¬ 
mined  at  each  E-layer  omtrol  point  from 

J  *  (1  +  0.0037  cos  (0.881  x)^-’  (15) 


where 


^12  -  12-month  running  mean  of  the  sunqiot  number 
X  s  sun’s  zenith  angle  -  degrees. 

The  value  of  the  absopdon  index  at  night  has  a  value  that  is  a  functiem  of  solar  activity  (Wakai,  1971) 
and  is  given  by 


Iain  “  0.025  (1  +  0.0311^2)  (16) 

Thevalueof I  usedisthemaximumofthevaluesgivaibyequations(15)and(16).Theabovevalue 
of  I  is  related  to  the  E(2(X)0)MUF  by; 

^(2000)«OT’  =  3,41  +  38,43*1  -  68.07*!*  89.97*1^ 

-  70. 97*I«  +  29.51*15  -  4.99*1*  (mHz)  ^  ^ 

The  f^E  is  determined  from  E(20(X))MUF  by  multplying  it  by  die  E-layer  distance  faeUx  for  zero 
distance  (0.208S).  The  minimum  of  the  values  at  the  E-layer  control  points  (^o^n  )  is  used  to  deter¬ 
mine  the  E-layer  MUF;  finally,  the  E-layer  MUF  (EMUF)  is  found  from 

EMUF  =  f^in  *  sec  (18) 


where  sec  is  the  secant  of  the  angle  of  incidence  of  a  ray  path  with  E-region  ionization  hei^t 

(110  km)  fer  E-modes. 
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The  minimum  hop  jnedicted  for  the  E-Layer  is  calculated  for  a  peak  layer  height  of  110  km 
and  die  padi  lengtih.  Tte  radiaticm  angle  associated  with  this  mode  must  have  an  angle  greater 
dum  some  prectetermined  minimum,  not  dictated  by  die  user,  the  angle  is  assumed  to  be  at  the 
horizon  (i.e.,  zeco  degrees).  TVansmissimi  losses  are  calculated  fm*  B-modes  fw  frequencies  up  to 
(1.4*EMUF). 

The  minimum  F2-layo’  hop  must  be  above  s<»ne  i»'eset  value  for  the  take-off  angle  as  fcr  the 
E-layer  and  must  also  penetrate  the  E-layer  at  die  P2-layer  angle  calculated  by  using  the  F2  peak 
layer  height  as  calculated  in  HFTDA.  The  P2-iaytf  h(^  are  increased  until  the  mode  take-off 
angle  sadsEes  the  above  restrictimis.  The  E*layer  penetratitm  frequency  is  (^o^n  * 

1.05),  where  sec  is  dw  secant  of  the  angle  of  incidence  of  a  ray  path  with  D-region 
ionization  height  (100  km)  for  P-modes. 

The  mixed  mode  is  calculated  by  using  the  E-layer  height  (110  km)  as  the  peak  laytf  height 
to  determine  a  take-off  angle.  The  frequency  must  be  suj^xxled  by  die  E-layer  at  one  eiul  of  the 
path  and  penetrate  at  the  other  end  for  this  mode  to  be  possible.  One  E-layer  hop  is  permitted, 
and  dien  the  path  is  completed  with  F2-layer  hops.  The  number  of  F2-layer  hops  depends  <m  the 
padileogth. 

Only  (me  minimum  hop  is  diosen  for  each  mode  of  {vopagadcm  for  each  layer,  since  the 
antenna  patterns  are  <mly  chosen  by  operating  frequency  and  padi  length,  thereby  making 
discrimination,  between  a  1-hop  and  2-hop  path,  impractical  until  more  realistic  antenna  patterns 
are  used. 

The  basic  transmission  loss  for  each  mcxle  is  given  by 

Lt,  =  *  Li  *  *  Lg  *  Lt,  idB)  (19) 

where  the  basic  free-space  loss  is  given  by 

=  32.45  +  20  1ogf  +20  logP'  (dB)  (20) 


and 


f  s  operating  frequency  (MHz) 

P'  =  virtual  slant  range  (km) 

Li  s  ncm-deviadve  ionospheric  absorption  below  the  MUF  (dB) 

Lc  s  loss  correction  facttx:  for  E-mcxles  (dB) 
s  over-the-MUF  loss  (dB) 

s  ground-reflecticm  losses  at  intermediate  reflection  points  (dB) 

(takmi  as  2  dB  per  ground  reflection) 

-  excess  system  loss  to  allow  for  aurcxal,  sporadic-E  obscniration, 
and  other  losses  not  explicitly  included  in  the  predictions  (dB).  In 
HFTDA,  it  has  a  fixed  value  of  7  dB. 

The  n(m-deviative  icmospheric  absoption  equation  is  the  one  used  in  lONCAP  with  the  near 
specular  reflecticm  losses  calculated  for  E-layo'  modes  at  low  frequencies  (Headri(tc  et  al.,  1971; 
Ibters  et  aL,  1983).  The  n(m-deviative  ionospheric  abscxpticm  is  given  by 
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(677 .2)  nl  sec 
(f  +  +  10.2 


idB) 


(21) 


where 


^  s  the  average  absorpticm  index  taken  over  all  the  control  p<wts 
»  gryofirequency  at  1(X)  km 
A  -  number  of  hops. 

The  ntm-deviative  absorption  given  by  equation  (21)  is  an  average  value  for  F2-laya’  modes.  The 
effects  of  E-region  electron  density  non-deviative  absorption,  collision  frequency  were  averaged  in 
the  curve  fitting  process;  therefore,  a  loss  cmrection  factor  for  E-layer  modes,  at  low  frequencies 
with  near  specular  reflection  valid  for  frequencies  above  2.0  MHz,  is  calculated  by 

L.  -  1.359  *  8.686-log(-g£^j  (effl)  (22) 


The  equations  for  E-mode  and  F-mode  losses  assume  that  the  mode  goes  throu^  the  absorbing 
regitm  (true  height  of  reflection  above  95  to  100  km). 

The  over-the-MUF  losses  are  calculated  by  using  the  I%illips  (I%illi{»,  1958;  Wheeler,  1966) 
method  for  values  of  tiie  frequency/MUF  ratio  up  to  approximately  1.4  to  1 J,  (Spending  on  the 
ground  distance.  For  values  of  this  ratio,  greater  is  considered  in  tte  scatter  regim  (Joint 
Technical  Advisory  Committee,  1960).  The  developer,  Don  Lucas,  of  this  model  was  an  original 
developer  of  lONCAP  in  the  mid-1980s  and  mule  significant  improvements  in  the  over-the- 
MUF  loss  model  when  he  develq)ed  the  code  for  HFTDA  in  1990.  While  HFTDA  uses  several 
models  from  lONCAP,  the  ovo'-the-MUF  losses  are  a  significant  dq>arture.  The  over-the-MUF 
losses  are  a  function  of  season,  sunspot,  geomagnetic  latitude,  local  time,  and  path  length.  The 
losses  are  calculated  as  follows;  first,  an  initial  ^ue  of  Lm  is  found  from  equation  [23];  where 

is  the  mode  MUF;  o  is  (0.0391  fm )  fix'  E-la^  modes,  and  (0.1172  )  is  for  F-layer  modes 

and  mixed  modes,  then 


L„  =  10logP  idB) 

o 


H  the  frequency  is  above  ,  the  frequences  at  which  scatter  effects  apply,  then  the  losses  are 
calculated  from 
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L,  =  10  logP +7.5  1og|-p^j  (dS) 

M 

P  =  — i—  f  exp/--^x*\  dx 
/Zni  \  2  / 

a 

The  parameter  i^and  is  found  from 

^•nd  “  for  d  s  1000  /cm 

^.Bd  =  (l-4  +  0.0175(-^'^-ggjj  *  f.  for  d  >  1000  /cm 

The  parameters  f«  and  0  are  defined  as  above. 

The  resulting  field  strength  for  each  mode  in  toms  of  the  system  loss  equatimi  is 

Ef  *  107 .2  +  Pj  +  20  logf  -G,  -L,  (dB  above  1  pv/m) 
Ef  =  107 .2  +  Pt  +  Gg  +  20  logf  - 


(24) 


(25) 


(26) 


with  the  system  loss  given  by 

L,  »  -  (Gf.  +  G^)  (dS) 


(27) 


where 

Pf.  a  effective  radiated  power  in  same  units  as  received  power  (Watts) 

bg  s  system  loss  (dB) 

G'c  =  antenna  power  gain  relative  to  isotropic  in  free  space  (dB)  for  the 
transmitter  antenna 

Gj  s  antmina  power  gain  relative  to  isotropic  in  free  space  (dB)  for  the 
receive  antenna 

The  systmn  loss  of  a  radio  circuit  is  defined  as  the  signal  power  in  decibels  that  is  available 
at  the  receiving  antenna  terminals.  This  excludes  any  transmitting  or  receiving  antenna  transmis¬ 
sion  line  losses,  since  such  losses  are  ccmsidered  readily  measurd>le.  The  system  loss  does 
indu&  all  die  losses  in  the  transmitting  and  receiving  antmina  circuits — not  only  the  transmis- 
si<m  loss  caused  by  radiation  from  the  transmitting  antmina  and  re-radiation  from  the  receiving 
antenna,  but  also  any  ground  losses,  dielectric  losses,  antenna  loadiag  coil  losses,  and  terminat¬ 
ing  resistor  losses.  Antenna  gain  is  taken  as  antenna  power  gain  that  is  the  product  of  antenna 
directive  gain,  in  the  direcdon  aligned  with  the  propagation  path  in  both  elevation  and  azimuth, 
and  of  antenna  effidency. 

Antenna  gains,  and  ,  are  in  the  direction  of  die  propagation  path  and  include  all 
antenna  losses,  sudi,  that  ^c4^r  is  an  ai^oximation  of  die  gain  Gp.  The  values  of  and 
are  required  for  any  elevation  angle,  azimuth  direction,  and  frequency.  In  HFTDA,  the  antenna 
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gains  are  determined  from  a  table  of  antenna  gains  for  each  type  antenna  represented,  which  is  a 
function  of  frequency  and  elevation  angle  (21  frequencies  and  6  elevation  angles).  After  the 
take-ofr  angle  is  determined  for  the  mode,  the  nearest  three  elevation  angles  in  the  table  are  used 
to  select  table  antenna  gains  for  the  nearest  integer  frequency,  above,  and  below,  the  desired 
frequency.  Then  three-point  Lagrangian  interpolation  is  used  to  interpolate  the  gain  for  the  two 
integer  frequencies.  The  results  of  the  two  Lagrangian  interpolations  are  then  interpolated  by 
using  a  linear  interpolation  based  on  the  given  frequency  at  the  two  closest  integer  table 
frequencies. 

After  the  field  strength  is  determined  for  each  of  the  three  modes,  the  predicted  Held  strength 
is  determined  from  the  sum  of  the  fields  in  microvolts  per  meter.  The  resulting  field  strength  is 
converted  back  to  dB  relative  to  1 -microvolt  per  meter. 

HFTDA  program  is  highly  modular,  allowing  changes  in  the  predicted  parameter  with  some 
ease. 

HFTDA,  called  LTLFLD  in  NOSC  Technical  Document  1848  (Systems  Exploration  Inc., 
1990),  has  been  developed  around  analytical  methods  and  empirical  data  for  the  HF  spectrum 
between  3  to  30  Mhz.  Any  calculations  or  predictions  outside  of  this  range  should  be  suspect. 

Sunspot  numbers  (SSNs)  over  150,  a  very  unlikely  event  for  the  next  20  years,  should  be 
used  with  caution  since  the  critical  frequencies  of  the  F2-layer  )  are  more  related  to  where 
you  are  within  a  given  cycle,  radier  thu  the  absolute  sunspot  number.  Some  scientists  believe  a 
saturation  effect  occurs  at  SSNs  of  greater  than  ISO  because  the  foF2  ceases  to  increase.  The 
impact  of  saturation  with  SSN  on  the  determination  of  transmission  loss  and  field  strength 
itself  is  not  clear. 

HFTDA  predictions  are  hourly  medians  of  the  monthly  medians  correlated  with  monthly 
medians  correlated  with  a  monthly  median  running  average  sunspot  munber.  Much  care  should 
be  taken  when  predicting  for  a  given  day,  or  a  few  days  ahead.  Field  strength  predictions  are 
weakly  associated  with  daily  values  of  solar  activity. 

Solar  flux  data,  that  are  used  as  an  indicator  of  ionization  in  the  F2-layer  may,  or  may  not, 
effect  in  the  same  manner  the  E-layer  and  D-lay^,  which  controls  most  non-deviative  absorption 
calculated  by  this  method.  Correlation  coefficients  are  not  available  for  the  D-,  E-,  and  F2-layers 
to  warrant  daily  predictions  of  field  strength.  The  accuracy  of  the  p>rediction  routine,  therefore, 
depends  entirely  upon  its  intended  use,  which  includes  frequency  assignment  in  the  long  term, 
siting,  antenna  selection,  day-to-day  frequency  use,  and  absolute  signal  determination.  The 
method  for  predicting  field  strength  at  a  receiving  location,  using  HFTDA,  is  the  solution  of 
Norton’s  transmission  loss  equation  (Norton,  1959).  It  is  this  equation  that  is  the  genesis  of  most 
analytical  approaches  to  HF  field  strength  predictions. 

The  simplest  calculations  were  used  to  increase  speed  of  the  calculation,  as  accuracy  was 
assumed  not  to  suffer  for  median  values  and  field  strength  calculations  alone.  Fourier  series 
expansions  representing  the  predicted  ionospheric  coefficients  (e.g.,  foF2 )  are  not  included  in 
any  of  the  variables. 

THE  lONCAP  MODEL 

The  lONCAP  HF  prediction  model  has  had  a  checkered  and  controversial  history  over 
its  first  decade.  When  first  delivered  to  the  U.S.  Army  in  1978,  it  was,  to  a  large  extent. 
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undocumented.  Ova*  the  ensuing  years  bugs  became  obvious,  with  several  analysts  working 
independently  toward  their  elimination.  Toward  the  late  1980s.  there  were  sevo-al  vmions  of 
lONCAP,  ea^  with  diffo’ent  modifications  and  improvements.  Attempts  wo-e  made  to  get  all  of 
the  changes  into  cme  version.  A  version,  called  VOACAP,  was  develc^ed  by  the  Voice  of 
Amoica  (VOA)  and  was  released  in  ^ril,  1993  (Lane,  Rhoads,  &  D^lasio,  1993;  Sweeney  et 
al.,  1993).  VOACAP  will  be  the  object  of  the  PENEX  field  strength  testing,  as  well  as  the  latest 
versicm  of  lONCAP  from  the  Institute  for  Telecommunications  Science.  The  following  descry)- 
tion  of  the  field  strength  calculation  is  taken  from  a  1983  lONCAP  Users  Manual  (Teters  et  al., 
1983). 

In  lONCAP,  die  basic  parameter  models  that  fmtn  the  basis  fOT  the  ionospheric  predictions 
are  numerical  map  r^resentations  of  the  parameters  that  describe  the  parameter,  their  temporal 
and  global  variations,  in  terms  of  Fourier  harmonics.  The  critical  fipequency  of  the  E-layer.  f^E , 
is  due  to  Leftin  (1976).  The  representation  of  the  critical  frequency  of  the  sp<M:adic-E  layer, 
fj^s,  is  due  to  Leftin  et  al.  ( 1968).  The  D-region  ionization  is  included  as  an  exponential 
decrease  bdow  the  E-layer;  a  height  of  maximum.  hmE,  of  110  km  and  hmEfYmE  ratio  of  S.S  are 
assumed  (i.e.,  Y^E  =  20  km).  lONCAP  uses  the  numerical  mtq)  of foFl  produced  by  Rosich  and 
Jones  (1973)  to  predict  foFl.  For  the  Fl-layer,  the  height  of  maximum  is  linearly  rdated  to  the 
sun's  zenith  angle,  and  the  ratio  of  the  height  of  maximum  hmFl  to  the  semi-thickness  of  the 
Fl-layer  is  assumed  to  be  4.  Models  of  the  critical  frequmicy  of  the  F2-layer,  foF2 ,  and  the 
monthly  median  M(30(X))F2  factw  are  due  to  Jones  et  al.  (1969).  The  height  of  the  maximum 
ionization  in  the  F2-region,  /4F2,  is  determined  from  the  M(3000)F2  factor  (Shimazaki,  1955). 
Having  obtained  hmF2,  die  semi-thickness  of  the  F2-layer,  YmF2,  is  obtained  frmn  4nF2  in  units 
of  YmF2  mapped  as  a  function  of  geomagnetic  latitude  and  solar  activity  for  solar  activity  (Lucas 
&  Haydon,  1966).  Numerical  maps  rqiresent  the  continents  for  use  in  ground-loss  calculations 
(Zadbiarisen,  1972). 

The  ionosphoic  profile  is  divided  into  four  regions:  D-E  region,  F2-region,  E-F2  valley,  and 
an  FI  ledge.  The  nose  of  the  E-layer  is  parabolic.  The  ionization  is  assumed  to  decrease 
exponentially,  starting  at  the  lower  part  of  the  E-region  (i.e.,  h,nE  -  0.85  ^E),  with  constants 
chosen,  such  that  the  slope  of  the  profile  is  continuous  at  this  point.  The  I^-region  is  assumed 
parabolic.  In  the  E-F2  region,  only  the  total  density  is  modeled  and  it  is  represented  by  a  line 
from  a  frequency)^  (0.98  fJE )  to  a  frequency  fy  (0.8516  )  at  the  E-layer.  The  Fl-layer  is 

considered  to  be  a  ledge  from  &e  F2-layer  to  tte  E-F  valley.  The  FI  ledge  may  be  either  a  linear 
layer  or  a  parabolic  layer.  If  the  height  of  maximum  ionization  of  the  FI  ledge  is  less  than  the 
height  of  the  F2-layer  at  the  frequency  of  foFl,  the  parabolic  shape  is  used.  If  the  height  of 
maximum  of  the  FI  ledge  is  greater  than  the  height  of  the  F2-layer,  then  the  height  of  maximum 
of  the  FI  ledge  is  reduced  to  the  F2-layer  height. 

lONCAP  uses  a  simple  closed  form  method  to  find  ray  paths  from  one  electron  density 
fx’ofile.  These  ray  paths  are  described  by  the  operating  frequency  f,  the  take-off  angle,  the 
vertical  height  of  reflection  h',  the  true  height  of  refiection  h,  and  the  ground  range. 

For  lONCAP,  the  basic  transmission  loss  for  each  mode  is  given  by 


*  Li  *  *  L„  *  +  Lp  *  *  Lg  (dB)  (28) 
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where 


the  basic  free-space  transmission  loss  expected  between  ideal,  loss-free, 
isotropic,  transmitting  and  receiving  antennae  in  free  space  (dB) 
non-deviative  ionospheric  absorption  below  the  MUF  (dB) 
loss  ctxrection  factor  fm*  E-modes  (dB) 
over-the-MUF  loss  (dB) 

ground  reflection  losses  at  intermediate  reflection  points  (dB) 
excess  system  loss  to  allow  fcxr  auroral  obscuration  and  other  losses  not 
explicitly  included  in  the  predictions  (dB) 
deviative  abscxption  (dB) 
sporadic-E  obscuration  loss  (dB) 
sp<H«dic-B  reflection  loss  (dB) 

The  basic  free-transmission  loss  ^bf  is  as  given  in  equation  (20).  The  non-deviative  ionospheric 
absorption  is  given  by  equation  (21)  with  the  absorption  index  at  each  control  point  given  by 

J  =  -0.04  +  exp  (-2. 937  +  0.8445  (29) 

and  the  gyrofrequency  fn  replaced  by  the  longitudinal  component  of  the  gyrofirequency  f  j  .  The 
formula  foe  absorption  index  J  is  in  terms  of  the  critical  frequency  of  the  E-layer  f^E ,  which 
includes  the  variation  in  zenith  angle  and  solar  activity.  This  formula  is  an  inversion  of  that  f(»merly 
used  frx  obtaining  fJE  from  I .  The  ctxrection  term  to  the  absmption  loss  for  E-layer  modes  that 
do  not  penetrate  the  whole  of  the  absorbing  layer  allows  for  frequencies  below  2  MHz.  It  is  similar 
to  equation  [22]  and  is  given  by 


=  A  +  B  logX^  (dB) 


(30) 


with 


A  =  1.359  fox  f„E  >  2  MHz, 

(  f„E  -  O.S]  ^ 

=  1.3591  °  ^  ^ - J  fox  0.5  s  fpF  S  2  MHz, 

=  0.0  fox  f„E  <0.5  MHz, 

B  =  8.617  fox  F^E  >  2  MHz, 


8.617 


(^) 


fox  0.5  s  fgE  S  2  MHz, 


0 . 0  for  f„E  <0.5  MHz, 


fy/  fgE  fox  h  >  90  km, 
fy  (90)  /  fgE  for  fr  s  90  km. 


true  height  of  reflection, 
equivalent  vertical  sounding  frequency. 


(31) 
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Tlie  equations  for  the  E-  and  F-laya*  mode  losses  given  by  equations  (29)  and  (30)  assume  that  the 
mode  goes  through  the  absorbing  region  (true  height  of  reflection  above  95  to  100  km).  When  the 
true  height  of  reflection  is  below  90  km.  these  equations  give  losses  much  higher  than  those 
observed.  In  lONCAP,  the  constant  term,  10.2.  in  equation  (21),  is  replaced  by  (V /2ir )  fcM*  values 
of  the  true  height  h  less  than  88  km  and  is  given  by 


=  63.07  exp  -2 


=  61+3 


I  h-7  0 
\  18 


^v-60) 

4.39  j 


(32) 


Because  a  complete  electron  density  is  used  in  lONCAP,  any  high-or  low-angle  mode  will  be 
considered.  Deviative  losses  are  considered  to  be  averaged  into  the  above  equations  for  reflec¬ 
tion  heights  less  than  that  at  the  layer  MUF.  For  modes  with  reflection  heights  greater  than  that 
at  the  layer  MUF  and  for  modes  just  past  the  E-P  cu^,  a  deviative  loss  term  is  added.  Tlie 
equation  for  this  term  Lp  is  based  on  the  relationship  that  the  loss  is  proptxtional  to  the  product 
of  collision  frequency  with  the  difference  between  group  path  and  phase  path.  The  deviative 
absorption  loss  in  dB  I>|,  is  given  by 


L  .  B(f^)  ih'  -  h)  „sec  (ipj  [{f^  +  +  10.2] 

^  *  (f  +  +  10.2 

where 

fv  = 
h'  = 
h  = 
f  = 
fi  = 

N  = 

= 

The  calculation  of  the  function  B(f  y )  dq)ends  on  the  presence  of  the  Fl-layer.  In  the  case  when 
the  Fl-layer  is  not  present,  fy  is  taken  as  fgE  for  the  determination  of  Lp  for  the  F2-layer,  to  pre¬ 
serve  continuity  at  the  E-  to  I^-laya*  transition.  In  the  case  when  the  Fl-layer  is  present,  fy  is  taken 
as  for  the  Fl-layer  deviative  abstnption  calculation,  and  fy  is  taken  as  foFl  for  the  F2-layer 

deviative  absorption  calculation.  This  assures  a  smooth  calculation  of  loss  for  all  electron  density 
profrles.  The  deviative  absorption  is  taken  as  the  maximum  value  of  all  the  possible  layer  hops  pos¬ 
sible. 

lONCAP  does  not  assume  that  propagation  ceases  at  the  MUF,  but  that  it  is  allowed  to 
decrease  acctxding  to  the  probability  that  the  critical  frequencies  of  the  E-,  F1-,  or  F2-layers  are 
above  the  predicted  values  of  the  MUFs  (Phillips,  1958;  Wheeler,  1966).  The  equation  for  this 
loss  tom  is 


vertical  sounding  frequency 

virtual  height  of  reflection 

true  height  of  reflection 

oblique  sounding  frequency 

longitudinal  component  of  the  gyrofrequency 

electron  density  profile 

angle  of  earth’s  ncnmal-to-ray  path  at  100  km 
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L,  =  10logi>  (dB) 

m 

(34) 

a 


where  f  is  the  oblique  sounding  frequency,  f ^  is  the  MUF  frar  the  circuit  elevation  angle  and  dis¬ 
tance,  and  o  is  the  Standard  deviation  of  the  distribution  of  the  MUFs  for  mode  being  calculated. 
Fw  the  E-  and  Fl-layers,  a  is  given  by  the  maximum  of  0.01  or  0.1  times  the  layer  MUF.  Fm*  the 
F2-laycr,  o  is  given  by 


1.28 


(35) 


whereFy  jiseitherFi  for/ :SF2MUF<»‘Fuf(X'y^F2MUF  from  table  2.  Table  2  shows  the  distribu¬ 
tion  of  daily  values  of  the  F2  MUF  about  their  monthly  median,  F2MUF,  as  the  ratios  of  upper  and 
lower  decile  MUFs  to  median  F2MUF  (f^  and  Fi ,  respectively)  fw  a  given  season,  a  given  solar 
activity,  4-hour  local  time  blocks  at  the  path  mi(4)oint,  and  for  each  10  degrees  of  geographic  latitude 
frmn  15  to  75  degrees,  north  w  south.  For  the  sporadic-E  layer,  o  is  given  below. 

The  ground  reflection  loss  calculation  f<x  randomly  polarized  sky  waves  assumes  equal 
amounts  of  energy  in  the  horizontally  and  vertically  polarized  fields.  The  losses  are  rq>resented 
by  the  following  equation 


=  10  lo,(.^4^)  (dB) 


(36) 


where  ICy  is  the  magnitude  of  the  vertical  reflection  coefficient  at  the  take-off  angle  A,  and  is 
the  magnitude  of  the  horizontal  reflection  coefficient  at  A.  These  r^ection  coefficients  are  a  func¬ 
tion  of  the  relative  dielectric  constant  of  earth,  the  conductivity  of  earth  (mhos/meter),  the  dielectric 
constant  of  free  space  (farad/meter),  the  frequency  transmitted,  and  the  take-r^  angle.  This  equation 
is  inadequate  for  take-off  angles  less  than  2  degrees. 

There  are  two  transmission-loss  terms  associated  with  modes  of  propagation  with  tiie 
sp<xndic-E  layer.  The  layer  is  modeled  as  a  thin  layer  occurring  at  the  height  h'Es  (usually  100 
to  110  km).  Its  effect  on  modes  of  propagation  passing  through  it  is  given  by  the  Es  obscuration 
loss.  It  is  calculated  by  a  method  proposed  by  Phillips  (1963)  and  modified  to  use  the  now 
available  mi^s  of  s,  the  median  value  of  the  highest  ordinary  wave  frequency  reflected 
from  the  spm'adic-E  layer  at  vertical  incidence: 
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-  10  logd  -  P)  (dB) 

z.LiIjEL 

a 

=  f^s  sec  4 


(37) 


The  variance  o  is  obtained  from  the  50%  and  10%  (or  90%)  values  of  f^s  acctxding  to 


f^s  f^s  -  SC 
1.28155  f^s 


(38) 


where  SC  is  the  tq)pa:  decile  of  f  s,  if  f  ^  fa  Es,  or  SC  is  the  lower  decile  of  f  o^s  iff  <  Es. 
PCX’  modes  that  have  reflected  firmn  the  sporadk-E  layer,  die  basic  transmission  loss  is  the  absorption 
losses^  c^^upplemented  by  a  reflection  loss  (corresponding  to  the  over-die-MUF  loss)  defined  by 

Lj,  =  8.91  P-®*’  (dB)  (39) 


Note  that  this  is  effectively  the  same  as  Lj,  ^  8.91  - 10  log  P .  The  probability  P  is  obtained  as  in 
the  above  equation. 

The  excess  system  loss  !<*  allows  for  auroral,  qxjradic-E  obscuration,  over-the-MUF  losses, 
and  other  losses  not  eiqilicitly  included  in  the  predictions.  In  lONCAP,  this  loss  is  determined 
from  table  3  for  path  lengtiis  less  than  2500  km,  and  table  4  for  padi  lengths  equal,  os  greater, 
than  2500  km.  llie  value  obtained  from  the  respective  table  is  the  median  value  for  a  given 
geomagnetic  latitude,  local  time,  and  seastm.  C^er-the-MUF  losses  are  subtracted  frrmi  the  table 
value.  When  the  sporadic-E  layer  is  calculated,  die  sporadic-E  obscuraticm  loss  is  subtracted 
from  that  value. 

The  sky-wave  field  strengdi  for  each  mode  is  directly  related  to  the  basic  transmission  . 
This  is  the  loss  as  contrasted  to  system  loss  diat  would  be  observed  if  the  actual  antennas  wav 
rqilaced  by  ideal,  loss-free  isotropic  transmitting  and  receiving  antennas.  The  field  strength  is 


E  =  107.2  +  20  log  fob  *  *  P  -I**  {isotropic)  (40) 


where 


E 


P 


rms  field  strength  in  dB  above  1  pV/m; 

transmitting  antenna  gain  (dB)  in  the  direction  of  the  ray  path  used 
to  determine  Lg  (decibels  referred  to  an  isotropic  anteima); 
transmitter  powo*  delivered  to  the  transmitter  antenna  in  decibels 
referred  to  one  watt; 
operating  frequency  in  MHz. 
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Table  3.  Expected  excess  system  loss  (dB)  for  padis  less  than  2500  km. 
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Table  4.  Expected  excess  system  loss  (dB)  for  paths  greater  dum  2500  km 


The  antenna  gains  in  lONCAP  are  determined  by  a  two-dimensicmal  liiuar  interpolatimi  <m 
a  table  of  antenna  gains  for  each  type  antenna  rqjresented.  The  table  gives  antenna  gains. 

46  elevation  angles  from  0  to  90  degrees,  for  every  integer  frequency  frcnn  2  to  30  MHz.  The 
resulting  array  is  46  by  29.  The  gain  determined  for  the  transmitting  antenna  im;ludes  a  factor  to 
accmmt  for  its  transmissicm  efficiency. 

The  fxedicted  field  strength  is  detomined  from  the  sum  of  the  field  strengths  of  die  individ¬ 
ual  modM  in  microvcdt  per  meter.  This  result  is  then  cmiverted  back  to  dB,  relative  to  1  micro¬ 
volt  per  meter. 

THE  ASAPS  MODEL 

The  ASAPS  model  (Ctfuana,  1993)  is  used  by  the  Australian  Government  Ionospheric 
Predicticm  Smvice  (IPS)  Radio  and  Space  Services  to  predict  HP  sky-wave  sy^em  p^ormance 
and  to  analyze  ionospheric  parameters.  ASAPS  stands  for  Advanced  Stand- Altme  Prediction 
System.  ASAPS  makes  good  use  of  gnqihical  displays  that  make  life  easier  for  the  user.  The 
model  inedicts  die  maximum  usable  frequency  (Ml^,  absoiptitm  limiting  frequency  (ALF), 
elevation  angles,  mode  (xobability,  padi  loss,  field  strength,  radio  noise,  noise-padi  locs(NP). 
and  signal-to-noise  power  ratio  for  a  range  of  possible  propagation  modes.  ASAPS  assumes 
great-drde  propagation,  widi  mirror  reflecti<»s  from  the  regular  E-  and  F-layers. 

In  ASAPS,  the  MUF  and  ALF  of  a  mode  are  used  to  (tetermine  the  availability  of  a  mode. 
The  MUF  is  evaluated  by  using  the  IPS  world  maps  of  ionospheric  characteristics  (Fox,  1988; 
Fox  &  McNamara,  1986, 1^).  The  mode’s  ALF  is  determined,  from  an  empirical  formula  in 
terms  of  the  solar  zenith  angle. 

The  computation  of  path  loss  is  based  oa  path  range.  For  drcuit  ran^  less  than  9000  km, 
the  basic  transmission  loss  is  essentially  the  method  described  in  the  Siq^lement  to  CCIR  Rqiort 
252-2  (COR,  1980).  For  paths  greater  dum  11,000  km,  the  path  loss  is  given  by  the  FTZ 
method.  For  path  Imigths  between  9000  aikl  11,000  km,  a  sq>arate  path  loss  is  computed  by 
using  bodi  tte  CCIR  mediod  and  the  FIZ  method.  The  resultant  path  loss  is  the  wetted 
average  of  these  two.  For  the  CQR  sq)arate-ioss-tam  procedure,  the  stq>s  in  the  computation 
indlude  (1)  determination  of  the  ray  pt^  of  the  mostly  normal  propagatitm  modes  between 
transmitter  and  receiver  on  purely  geometrical  terms,  (2)  piedicdon  of  the  state  of  the  iemosphere 
at  each  sky  reflection  point  to  determine  which  of  these  modes  are  open — the  q)erating  frre- 
quency  must  fall  between  the  ALF  and  the  iqq}er  decile  of  the  MUFs  frxr  the  path — and  (3)  com¬ 
putation  of  die  path  loss  for  eadi  open  mode.  The  FTZ  method  makes  circuit  predictions  fix 
only  a  composite  mode  (i.e.,  no  specific  modes  are  ermsidered). 

The  sqiarate  transmission  loss  and  gain  factors  that  are  specifically  taken  into  accoimt  in 
ASAPS  for  paths  lengths  <  11,000  km  are  as  follows: 

•  spatial  attenuation,  including  focusing  of  rays  with  low  elevation  angles  and  rays 
propagated  to  very  long  distances  and  G^) 

•  non-deviadve  and  deviadve  absorption  (I>^ ) 

•  polarization  coupling  loss  (Lp ) 

•  ground-reflection  loss  (L^ ) 

•  sp(xadic-E  obscuration  loss  (L^ ) 

•  transmitting  and  reedving  antenna  gains  (Cf  and  ,  respectively) 
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The  btsic  transmission  loss  for  path  lengths  <  11,000  km  is  given  by 


*  Lp  *  Lg  *  L„  - 


<dB) 


(41) 


The  basic  free  space  attenuation  is  the  same  as  given  in  equatitm  (20)  for  HFTD  A.  Tlw  gixnind 

refto^ioD  loss  Lg  for  mult4>le>hop  modes  is  the  same  as  givmi  in  equatitm  (36)  for  lONCAP. 
Depending  on  circuit  path  lengdi,  die  horizon  focus  gain  is  calculated  using  one  of  the  fr^owing 
two  methods: 

For  drmiit  paths  that  are  less  dian  one-quarto'  of  the  Earth’s  circumference  in  length  (about 
10,(X)0  km),  the  hisiz<m  focus  gain  for  E-  and  F>modes  is  calculated  by  the  Bradley  (1970) 
method.  This  method  gives  as  a  function  of  the  mean  elevation  angle  of  the  upgoing  rays 
over  all  hops .  It  is  indqpendent  of  tiie  number  of  hops,  and  it  reaches  a  maximum  of  9  dB  at 
grazing  incidence.  G,  has  been  approximated  by  four  6-degree  Chebychev  polynomials: 

•  For  mean  elevatitm  angles  ^  10  degrees 


Gf  {E-modes) 
GfiF-modes) 


»  4.14  -4.00x  +  3.007x^  +0.069x^ 

-  1 . 06x*  +  0 .681x5  -  0 . 347X* 

»  6.03  -  3. 161X  +  1.4X*  ■*•0.624x5 

-  1 . 413X*  ■*•  0 . 088x5  +  0 . 533x« 


(42) 


where 


X  =  Jtean  elevation  angle  -  5 


(43) 


•  Fa*  mean  elevation  angles  >  10  degrees 


Gf  (E-modes) 
Gf{F-modes) 


=  0.81  -0. 876x  0 . 353x5  _  q  _  028x5 

-  0 . 227x*  -  0 . 274x5  *  0 . 387 x* 

«  1.50 -1.244X■^  0.173x5  •*■0.540x5 

-  0 . 267x*  -  1 . 136x5  +  0 . 853x5 


(44) 


where 


_  mean  elevation  angle  -  50 
40 


(45) 


Fa:  drcuit  path  lengths  that  are  greater  than  aie-<iuarter  of  the  eartii’s  circumference,  the  horizon 
focus  gain  G^  is  given  by  the  following  formula  (Hotenbach  &  Rc^a,  1979): 


Gf  =  -20  lo 


4-^1) 


(46) 
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wlioe 


D  s  grett-drcle  path  length  between  transmitter  and  receiver  in 
kilometers 

^  B  lfarl>  <  diree-quarters  of  tlM  Earth's  circumference  (3a  R/2) 
and  2  for  longer  distances 
R  B  Eardi’s  radius  (6370  km). 

The  horizon  focus  gain  is  limited  to  a  maximum  value  of  30  dB. 

The  types  of  absoiptioo  taken  into  aocmmt  in  ASAPS  are  die  normal  nmi-deviadve  and 
deviadve  absorption  (including  the  median  winter-anomaly  absorption).  The  auroral  absorption, 
the  extra  attenuatkm  of  radio  waves  that  traverse  die  aurar^d  zones,  is  not  calculated  in  ASAPS. 
The  estimadmi  of  the  nmi-deviadve  and  deviadve  absorption  in  ASAPS  is  based  (» the  analysis 
of  vertical-incidence  absorption  measurements  (George,  1971).  The  absorptimi  estimatim  uses  a 
relationship  between  the  ateorptimi  at  vertical  and  oblique  incidence  obtained  by  George  and 
Bradley  (1973)  frmn  the  results  of  ray-tracing  calculations  throu^  model  ionoqiheres.  Tte 
procedure  is  described  in  detail  by  George  tad  Bradley  (1974).  This  takes  no  account  of 
absorption  in  die  F-region  arising  from  coUisuxu  between  electrons  and  ions. 

The  absorption  ,  eiqierienoed  by  the  ordinary  wave  frequency  ^ob  reflected  tdiliquely 
from  the  icmoqphere,  is  given  in  terms  of  the  absorption  (jf^ )  of  the  ordinary  wave  at  vertical 

incidence  on  a  related  frequency  f  ^ 


i-.(fob) 


X.^(fv)  ify  *  fjt)^seci^oo 
(^06  ♦ 


(47) 


where  f  j ,  the  electrrm  gyrofrequency  about  the  vertical  component  of  die  geomagnetic  field,  is 
taken  as  positive.  This  is  sfproximately  1.5  MHz  for  the  B>modes  and  1.0  MHz  for  die  F-modes. 
The  parameter  fy  is  given  in  toms  of^ojb  by 

fy^f^cosi^gg  (48) 

iloo  is  the  inddence  angle  of  the  unrefracted  ray  at  a  layer  hei^t  of  100  km,  such  that 

sini^gg  s  0.9845CO8A  (49) 

wtere  A  is  the  ray’s  elevation  angle  at  the  ground. 

George  (1971)  showed  that  at  no(m  the  term  A(fv)  bL^  (fy)(fy+fj  )2  was  related  to 
dw  limiting  value  of  A  ( f  ^1  for  a  sufficiently  hi^  frequency  that  signals  traverse  the  whole  of 
the  absorbing  tegimi  without  deviation,  and  is  a  function  of  the  ratio  of  Jfv  to  f^E  as  given 
by 


(50) 
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The  funcdon  is  i^^oximately  independent  of  locati(».  seastm  or  solar  epoch  and  is 
given  in  figure  1.  (The  figure  given  1^  ai^  used  by  ASAPS  is  the  same  as  figure  7  in  the 
Supplonmit  to  CCIR  Rqnxt  252-2  (CCIR,  1980),  exc^t  that  the  maximum  value  of  is 
liniited  to  1.56.)  Equation  (50)  holds  under  normal  absorption  conditimis  and  also  aj^lies 
ai^oximately  Ifor  winter-anomaly  absorpticMi.  Samuel  and  Bradley  (1975)  have  shown  dtat 
equatkm  (50)  may  also  be  extended  to  other  times  oi  day  than  norm.  For  a  given  m(»di  and 
location.  A,  increases  linearly  with  12-m(mth  running  mean  sunspot  number  ^12  and  changes 
diumally  as  a  function  of  the  solar  zenith  angle  X-  For  a  given  m<mth  and  location.  A,  in 
ASAPS  is  given  by 


Ar(T,x)  =  Aj.(0,0)  cos'(0.881x)  <1  +  0.0067  D  (51) 


where  A,.  (0,0)  is  die  value  of  A,  at  solar  minimum  (T  sO)  and  solar  no<m  (JL  s  0).  The  parameto’ 
T  isthei<»ioq>hericin(kixusedbyIPS,andit8relati(»ish4>tol?x3  isgivenintable5Cnimer,  1968). 


Hgure  1.  The  George  absorption  function  . 


Table  5.  RelaticMiship  between  the  IPS  ionosphoic  index  T  and  R12 
(Tiirner,  1968). 


Mcmth 

Ordinates  at  origin  (c) 

Slq)e  (d) 

January 

-1J8 

0.96 

FbtMTUI^ 

0.81 

0.94 

March 

-0.65 

0.97 

^ril 

-054 

0.99 

May 

-0.06 

1.00 

June 

0.34 

1.01 

July 

-0.18 

1.00 

August 

1.19 

0.97 

Sqjtember 

2.70 

0.93 

Ot^ber 

1.87 

0.93 

November 

2.45 

0.92 

December 

-0.04 

0.96 

where 


T  =  c  +  d  *  Ri2 

The  parameters  (0,0)  and  p  s's  functions  of  month  and  of  modified  dip  latitude  /  X  / 
where 


tanJT  » 


I 

v/cos  A 


(52) 


The  parameter  j  is  the  magnetic  dip  angle  in  radians  (considered  positive  if  north  of  the  magnetic 
equator) ,  and  A  is  the  geographic  latitude  of  the  point.  Figures  2  and  3  (figures  8  and  9  in  the  S  tq^le- 
ment  to  CCIR  Repot  252-2  [COR,  1980])  give  p  and  Aj.  (0,0),  respectively,  for  each  month  and 
/  X  /  value.  At  high  latitudes  with /  X  /  >70°,  values  of  A,  (0,0)  and  p  are  taken  as  for  /  X  /- 
70°;  hence,  combining  equations  (47),  (50),  and  (51)  gives  ^  daytime  absoption  for  ASAPS 


Ay(0,0)cos^(0.881x)  (1  +  0.0067 


"  M-6) 


sec  2 


100 


(53) 


For  a  multiple-hop  mode,  the  absoiption  on  each  hop  is  evaluated  sq)arately  by  using  mid-hop  posi¬ 
tion  icmospheric  characteristics  and  a  mean  elevatioi  angle  fo  the  upward  and  downward  legs  of 
the  hop. 

After  sunset,  the  absorption  falls  to  a  small  non-zero  value.  The  total  nighttime  absorption 
(Wakai,  1975)  is  given  by 


_  (7  +  0.019D)  (1  +  0.015D 


iO 


(ds) 


(54) 
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where is  the  (^)erating  frequency  in  MHz,  and  is  the  totr  1  ground  range  between  transmitter 
and  receivor  in  km.  The  values  determined  are  regarded  as  iq^lying  sq>arately  to  each  of  the  pn^a* 
gatitm  modes  that  can  exist.  The  value  fa  ,  detennined  from  equation  (54),  is  taken  for  the  mode 
when  it  gives  a  value  larger  than  that  obtained  from  equation  (S3). 

When  an  upgoing  wave  is  incident  on  the  ionosphere,  it  leads  to  the  excitation  of  an  (X’diruuy 
(O)  and  an  extraodinary  (X)  wave.  These  waves,  having  different  polarizations,  may  be 
r^arded  as  propagating  independently  within  the  ionosphere  and  being  subject  to  ddferent 
amounts  of  absoption.  This  polarization  coupling  lossL^  is  given  in  ASAPS  by  the  empirical 
formula 


Lp  =  A  +  fl*  (hop  length)  *  inumbez  of  hops)  (dB)  (55) 

where  the  cmistants  A  and  B  depend  on  the  mode,  as  defined  in  table  6.  and  the  h(^  length  is 
espessed  in  radians,  with  a  maximum  value  equivalent  to  a  hop  length  of  1S,(X)0  km. 

Table  6.  Polarization  loss  constants. 


Qmstants 

E-modes 

F-modes 

A 

2.400 

B 

2.548 

1.019  1 

Sporadic-E  obscuration  losses  are  taken  into  account  in  ASAPS  by  using  an  empirical 
frxmula  developed  by  Sinno  et  al.  (1976).  Pex*  a  wave  frequency /in  MHz  and  an  incidence 
angle  iioo  of  the  oblique  ray  at  a  height  of  110  km,  the  obscuration  loss  for  one  traverse  of  the 
Es-layer  by  an  F-mode  is  given  by 


L,  =  -log  (1 -S’*)  (dB) 


where 


S 


(f^sseci,,o] 


1  +  10 


(56) 


(57) 


The  sporadic-E  critical  frequency  f^E  s  is  approximated  by  1.2  *  ff,E ,  with  a  lower  limit  of 
1.2NQlz.Equation(56)givesL^=O.O4dBf<X'/s/0f  s^seciiocandL^  -*  »  as/-*0. Because 
sporadic-E  obscuration  loss  can  become  quite  large  for  small  frequencies,  sp^adic-E  obscuration 
loss  in  ASAPS  is  limited  to  a  maximum  value  of  54  dB  per  hop.  is  evaluated  separately  for  each 
hq},  then  summed. 

Finally,  for  circuit  lengths  <  90(X)  km,  ASAPS  determines  the  rms  field  strength  from 


Ej^  =  107 .2  +20  logf  +  (dB  above  1  pV/m)  (58) 
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Pt  =  signal  powo*  (in  dBW)  at  the  ii^ut  to  the  transmitting  antenna 
f  =  transmitting  frequency  in  MHz 

>  transmitting  antenna  gain  in  dB  relative  to  an  isotropic  antenna  in 
die  direction  of  the  prt^agation  mode 
s  basic  transmission  loss  fot  the  mode  in  dB 

The  antenna  gains  that  are  in  ASAPS  are  determined  by  a  two-dimensional  linear  interpolatitHi  frmn 
tables  of  antmma  gam  given  at  6  take-off  angles  and  21  frequencies. 

For  circuit  lengths  >11,CXX)  km,  ASAPS  uses  a  variance  of  the  FTZ  method  adopted  by  the 
CQR  for  use  at  circuit  ranges  greater  than  9000  km  (CCIR,  1990a).  The  median  field  strength  is 
given  by 

^m»dimn  “  -66.4  *  -  Ly  (dB  above  1  \i.v/m)  (59) 

where 

»  139.4  -  20  logU' 

X .  1  1  («» 

<  » ♦  ( fl.'-f .) »  ” 


fiee-space  field  strength  fm  3,000  kW  eip 
virtual  slant  path  length  in  km  between  transmitter  and  receiver 
upper  limit  frequency  (operational  MUF)  given  in  equations  [2] 
a^  [3]  with  1.2  replacing  1.0  MHz 

lower  limit  frequency  (MHz)  (similar  to  equation  [4],  but  includes 

winter  anomaly  effects 

F-layer  gyrofrequency  (MHz) 

transmission  frequency  (MHz) 

transmitter  power  in  dBW 

antenna  gain  of  transmitter  antenna  in  dBi 

focus  gain  on  very  long  paths  in  dB 

extra  empirical  loss  (or  gain),  currently  defined  as  -4.2  dB 


The  focus  gain  is  given  by 


G.P  =  10  log  - (dB) 


(61) 
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where  R  is  the  Earth’s  radius,  and  D  is  the  great-circle  distance  between  transmitter  and  receiver, 
in  km.  G,p  is  limited  to  a  maximum  value  of  IS  dB,  since  it  tends  towards  infmity  as  D  approaches 
a  multiple  of  tc  i? .  llie  rms  Eeld  strength  is  taken  to  be  1 .6  dB  above  the  median  value 

-  1-6  (62) 


By  equating  the  CCIR  and  FTZ  rms  field  strengths,  the  equivalent  FTZ  pathloss  becomes 
^FTZ  *  L*  =  167 .8  +  20  logf  -  E^X  -  G,p  (dB)  (63) 

For  circuit  path  lengths  between  9,(XX)  and  11,(X)0  km,  a  separate  pathloss  is  computed  in 
AS  APS  by  both  the  CQR  method  (l^)  given  in  equation  [41]  and  the  FTZ  method  (L^  given 
in  equation  [63].  The  resultant  pathloss  (L^^  )  is  the  weighted  average  of  the  two 

=  (1  -  w)  Li,  +  wL^  (dB)  (64) 


where  the  weight  w  is  given  in  terms  of  the  path  length  (d  )  by 

^  _  D  -  9000 
2000 


(65) 


The  rms  field  strength,  for  circuit  path  lengths  between  9,(X)0  and  11,000  km,  is  then  given  by  the 
CQR  formula  in  equation  [58]  with  L^,  replaced  by  L^y . 

THE  AMBCOM  HF  PREDICTION  MODEL 

The  AMBCOM  (Hattield  &  Smith,  1987)  is  a  fast,  versatile,  two-dimensional  (2D)  ray 
tracing  program  for  simulating  high-frequency  (HF)  ionospheric  sky-wave  propagation  and  the 
resulting  performance  of  various  HF  systems,  for  example,  communications,  broadcast,  geoloca- 
tion,  surveillance  and  over-the  horizon  radar  (OTHR)  systems.  Its  name  is  an  acronym  for 
Ambient  Communications,  by  analogy  with  the  Nuclear  Effects  on  Communications  (NUCOM) 
code,  from  which  it  is  largely  derived. 

The  basic  premise,  underlying  the  design  of  AMBCOM,  is  that  the  accuracy  of  the  propaga¬ 
tion  model  should  be  compatible  with  the  accuracy  available  in  the  ionospheric  model.  Since 
only  enough  ionospheric  profile  information  is  normally  available  to  portray  three-parameter 
layers,  little  purpose  is  served  by  employing  a  ray  path  calculation  that  presumes  more  detailed 
knowledge  d  the  profile.  On  the  other  hand,  strong  horizontal  gradients  in  critical  frequency  and 
height  are  iy>parent  even  in  maps  of  monthly  median  parameters.  Neglect  of  these  systematic 
variations  can  be  equally  misleading:  gross  eirors  in  range  estimation  or  in  the  computed 
elevation  angle  at  the  receiver,  or  both,  can  occur  under  certain  conditions.  AMBCOM  is 
designed  to  account  for  such  variations  to  the  frrst  order. 

The  distinguishing  feature  of  AMBCOM  is  the  prc^agation  model,  which  combines  the 
advantages  of  open-ended  ray  tracing  with  the  speed  of  analytic  integration  through  parabolic 
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layers.  By  “open-«nde<l”  it  is  meant  that  the  path  of  the  ray  through  the  ionosphere — and  the  range 
at  which  it  returns  to  earth — are  determined  by  the  refractive  effects  of  the  ionosphere.  This  is  in 
contrast  to  virtual-geometry  codes,  where  the  point  at  which  the  ray  returns  to  earth  is  specified  in 
advance  and  each  h(^  is  symmetric  about  its  predetermined  reflection :  .int.  Examples  of  widely 
used  virtual-geometry  codes  are  HFMUFES,  lONCAP,  and  RADARC.  Such  codes  cannot  predict 
‘‘unconventional”  modes  of  pr<^agation,  such  as,  tq>side  reflections  (M-modes)  and  chcx’dal  modes. 
Because  AMBCOM  uses  ray  tracing,  it  will  handle  large  ionospheric  gradients  and  will  predict 
asymmetric  h(^s  and  imconventional  modes  when  they  are  indicated  by  the  ionospheric  model.  Sig¬ 
nificant  differences  in  the  ionosphere  and  pr(^agation  models  used  by  AMBCOM,  RADARC,  and 
lONCAP  are  summarized  in  tables  7  and  8. 

Another  AMBCOM  feature  is  that  it  is  one  of  the  most  thoroughly  documented  prc^agation 
codes.  Although  originally  developed  for  mainfirame  computers,  a  PC-based  version  has  been 
develqjed  recently. 

AMBCOM  is  based  on  the  NUCOM  code  develqjed  by  the  Stanford  Research  Institute 
(SRI)  during  the  1960$  and  1970s  under  sponsorship  of  the  Defense  Atomic  Support  Agency 
(DASA)  and  its  successor,  the  Defense  Nuclear  Agency  (DNA),  to  predict  the  p^ormance  of 
HP  communication  systems  under  normal  (i.e.,  ambient)  and  nuclear  ionospheric  conditions 
(Nielson  et  al.,  1967).  AMBCOM  employs  the  ray  tracing  and  communication  system  concepts 
of  NUCOM,  but  it  is  intended  primarily  for  use  under  ambient  conditions.  The  current  version  of 
AMBCOM  incorporates  additional  mo^ls  of  the  ambient  ionosphere  that  were  developed  in  the 
late  1970s  (Hatfield,  1980). 

Table  7.  Summary  of  significant  model  diffo-ences. 


Models 

Program  Names 

RADARC 

AMBCOM 

lONCAP  1 

1  Ionosphere  Generation  | 

Median  models 

Yes 

Yes — including 

Auroral  Ionosphere 

Yes 

Spatial  representation 

4  samples 

41  samples 

4  samples 

Real-data  input 

ionograms  at  4  loca¬ 
tions 

9  parabolic  parame¬ 
ters  at  up  to  41  loca¬ 
tions 

ionograms  at  4  loca¬ 
tions 

|Prop{^tion  Model 

Raytradng  method 

Martyn’s  Theorem 

Semi-analytic  raytrace 

Martyn’s  Theorem 

Tilts,  gradients 

No 

Yes 

No 

Topside  reflections 

No 

Yes 

No 

Radar  propagation 

Yes 

Yes 

No 

Point-to-point  propa¬ 
gation 

No 

Yes 

Yes 

1  Radio  System  Simulation  | 

1  Circuit  reliability 

No 

No 

Yes  1 

34 


Table  8.  Differences  in  the  input  ionospheric  and  noise  data. 


Program  Name 

Variables 

RAD  ARC 

AMBCOM 

lONCAP  1 

m _ 1 

Coefficients 

reference 

Oslo  (CCIR,  1967) 

New  Delhi  (Jones  & 
Obitts,  1970) 

Oslo  (CCIR,  1967) 

-  SSN  variation 

Linear 

Second  order 

Linear 

-  Centered  on 

15th  of  month 

15th  of  month 

15th  of  month 

foFl 

Analytic  (formula) 

Filled  layer  based  on 
E,F2 

Coeffidents  (Rosich, 
1973) 

foEs 

Coefficients  (Leftin 
et  al.,  1968) 

Coefficients  (Leftin 
etal.,  1968) 

Coefficients  (Leftin 
et  al.,  1968) 

foE 

Coefficients  (Leftin, 
1967) 

Coefficients  (Leftin, 
1976) 

Coeffidents  (Leftin, 
1976) 

|M(3000)F2  1 

-  Method  used 

Coefficients 

Coefficients 

Coefficients 

-  SSN  variation 

Linear 

Linear 

Linear 

Land  mass 

Coefficients 

Coeffidents 

Coefficients 

[Atmospheric  Noise  | 

-  Coefficients 

Yes 

Yes 

Yes 

lime  variable 

LMT  (hours)  (Lucas 
&  Harper,  19^) 

UT  (hours)  (Zachari- 
sen  &  Jones,  1970) 

LMT  (hours)  (Lucas 
&  Hai^r,  1965) 

The  ionosphere  is  modeled  with  three  parabolic  layers  of  electron  density.  The  £^^2  is 
represented  by  a  set  of  coefficients  develq)ed  by  Jones  and  Obitts  (1970).  The  M(3000)F2  is 
represented  by  coefficients  due  to  Jones  et  al.  (19^).  The  height  of  the  F2-layer  is  determined 
from  the  M(3000)F2  by  using  a  relationship  by  Shimazaki  (1955).  This  result  is  corrected  fw 
retardation  in  the  E-  and  Fl-regions  (Wright  and  McDuffie  (1960).  The  is  reprwented  by 
coefficients  due  to  Leftin  et  al.  (1968).  The  height  of  the  E-layer  is  modeled  at  115  km,  with  a 
semi-thickness  of  25  km,  and  with  no  day/night  variation.  The  height  and  semi-thickness  of  the 
pseudo-Fl  layer  are  computed  fi-om  the  parameters  of  the  E-  and  F2-parabolic  layers.  The 
bottom  of  this  layer  is  set  at  130  km,  and  this  layo:  is  required  to  overlap  with  the  F2  lay^  by 
one-half  its  own  semi-lhickness  (75  km).  The  aitical  frequency  jf^Fl  is  related  to  by 

f^Fl  = 
n  - 

120 


The  maximum  value  of  is  limited  to  being  less  than  0.695  fj?2 .  The  f^s  is  represented  by 

coefficients  developed  by  Leftin  et  al.  (1968).  The  height  and  semi-thickness  of  the  Es-layer  are  set 
to  110  km  and  1  km,  respectively. 
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special  features  of  AMBCOM  include  the  following: 

•  A  model  of  the  electron  density  profile  in  the  high  latitude  ionosphffe  (including  the 
aurwal  z(x»). 

•  A  model  for  computing  auroral  absorptimi. 

•  Models  for  computing  reflectimi  ami  penetration  losses  fm  the  Es-layo',  as  fimcticms 
of  the  maximum  frequency  and  blanketing  frequency  of  the  layer. 

•  Cmisida^on  of  both  topside  and  bottomside  reflecticnis  from  die  Es-layer  (e.g.,  M 
mottes). 

•  An  option  to  search  for  high  rays  on  one-  and  two-hop  F2  modes  (IFIH  and  2F2H 
mod^). 

The  new  model  of  d;  nigh-ladtude  ionosphere  is,  itself,  a  modification  of  the  RA£>C- 
POLAR  model  developed  by  Elkins  and  Rush  (1973a,  1973b).  It  is  based  on  the  incoherent-scat¬ 
ter  radar  measurements  of  electron  density  profiles  in  the  auroral  ionosphere  (Vondrak  et  al., 
1978)  by  SRI  at  Chatanika,  Alaska.  The  auroral  absorption  model  was  devel(^>ed  by  SRI  from 
previously  compiled  riometer  data.  It  describes  the  variation  of  median  auroral  abs^tion  as  a 
function  ^  (xmected  geomagnetic  latitude,  longitude,  season,  local  time,  solar  activity,  and 
geomagnetic  activity  (Vondrak  et  aL,  1978). 

The  expression  for  basic  transmission  loss  foe  AMBCOM  is  given  by 


Lb  -  Lt,f  +  Ljjs  *  Lj^i  *  LJ^2  +  Lf^  +Lp  +  *Lg  (dB)  (67) 

where 

Lbf  =  free-space  loss  (dB) 

Ld2  =  divergence  loss  for  F2  high-ray  modes  (dB) 

Lju  =  nuclear  D-region  absorption  (dB) 

La2  s  E-Fl  absoiption  (dB) 

La3  =  non-deviative  D-region  abstxption  (Li)  (dB) 

Lp  =  F1-F2  deviative  absoption  (dB) 

Laux  =  aurmal  absoq)tion  (dB) 

Lq  =  SpOTadic-E  obscuration  loss  (dB) 

Lg  =  Spcx-adic-E  reflection  loss  (dB) 

Lg  =  summed  ground-reflection  loss  at  intermediate  ground-reflection 

points  (dB). 

The  free-space  spreading  loss  term  fm  AMBCOM  is  the  same  as  that  given  by  equation  (20).  The 
d^ocusing  loss  term  foe  F2-layer  high-ray  modes  is  given  by 


Lo3  =  10  log 


R  sin  -~ 
R 


sinA 

cosP 


(P') 


dD 


idB) 


(68) 
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where 


=  group  path  (km) 
fi  s  take-off  angle 

A  =  arrival  angle 

D  =  ground  range  (km) 

R  =  eardi  radius  (km). 

To  improve  the  accuracy,  Lg^  is  calculated  for  bodi  forward  and  reverse  rays,  then,  die  average  of 
these  two  values  is  used.  For  the  reverse  calculation,  angles  A  and  P  are  interchanged. 

AMBCOM  provides  models  fv  estimating  several  components  of  abs(»ption  loss,  acceding 
to  the  altitude  at  which  the  absmption  occurs  and  die  source  of  the  ionization  responsible  fcH*  the 
absoiption — solar  radiation,  nuclear  radiadmi,  or  auroral  precipitation.  Hiese  are  computed 
individually  (if  apprqpriate)  in  various  parts  al  the  program  and  are  cmnbined  at  the  end  of  each 
hop.  Up  going  and  down  going  legs  of  a  hop  are  treated  separately  to  take  into  account  spatial 
variations  in  electron  density  and  possible  changes  in  the  angle  of  incidence  of  the  ray. 

Solar-amtrolled  absmption  is  conqiuted  in  three  height  regimes:  the  D-regitm  (below  92 
km),  where  the  absorption  is  non-deviative;  the  E-Pl  region  (92  to  ~  ISO  km),  where  both 
deviative  and  non-deviative  absoption  may  occur;  the  F1-P2  region  (above  150  km),  where  the 
absoiption  is  treated  as  stricdy  deviative.  (If  desired,  the  nuclear  absorption  can  be  computed  for 
the  same  three  height  regimes;  however,  the  nudear  routines  are  bypassed,  unless  externally 
generated  electron-density  protiles  are  iiqiut.)  The  solar  D-region  absmption  ^a3  for  the 
ordinary  ray  is  evaluated  at  each  D-region  traversal  by 


^A3 


119(1.0  +  0.0037  SSW)  (cos0.881x)^'^sec^c 
(f  + 


fox  X  i  102* 
for  X  >  102* 


(69) 


where 

SSN  s  sunspot  number 

X  s  solar  zenitii  angle  at  location  of  D-region  crossing 
=  ray  angle  of  incidence  at  70  km 
f  =  operating  frequency  (MHz) 
f„  -  electron  gyrofrequency  at  100-km  altitude  (MHz). 

The  second  solar-controlled  absoiption  region  extends  from  92  km  to  about  150  km  and  considers 
total  ionospheric  abstxption  in  each  parabolic  layer  (both  deviative  and  non-deviative).  The  same 
expressirm  is  used  fix  ambient  and  nuclear  cases,  but  the  values  of  some  of  the  parameters  differ. 
The  term is  computed  by 


A.Aj 


^  AA_£  dS/ layex  txavezsed 

f2 _ Ay  dB/siat. 

-  py(2ya  -  y) 
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where 


Ay 

y 

y. 

V 


4>o 

f 

fn 

fc 


aiKl 


slab  thickness  (km) 

vertical  distance  of  the  slab  from  the  laytf  botton  (km) 
layer  semi-thickness  (km) 

collision  frequency  with  air  molecules  and  oxygen  atoms, 

apfxopriate  to  the  slab  altitude  (millions  of  collisions/s) 

inddence  angle  at  the  entrance  to  the  layer 

opoating  frequency  (MHz) 

electron  gyrofrequency  at  the  slab  height 

layer  critical  frequency  (MHz) 


P  = 


(71) 


At  altitudes  greater  than  about  150  km,  the  collision  frequency  is  generally  about  three  orders  of 
magnitude  smallo*  than  the  wave  frequency.  This  leads  to  a  relatively  simple  e]q>ressi(»i  for  the  total 
deviative  absmption  incurred  as  given  by 

ip  =  (P'  -P)  (dB>  (72) 

where  P'  and  P  are  the  group  and  phase  paths  for  that  portion  of  the  ray  path  above  ISO  km,  and 
v  s  103  s*l.  The  deviative  absoption  is  usually  small,  excq>t  at  near-vertical  incidence. 

The  auroral  abstxpdon  model,  is  based  on  a  method  proposed  by  Foppiano  (1975)  and 
adopted  (after  certain  modifications)  by  the  CGR  (1980).  The  Foppiano  (1975)  model  was 
developed  by  fitting  equations  to  published  riometer  data  from  27  stations  in  the  northern 
hemisphere  to  model  the  spatial  and  temporal  variations  of  auroral  absmption.  The  resulting 
equations  can  be  used  to  predict  the  monthly  oiw-way-vertical  absorption.  Am.  at  a  frequency  of 
30  MHz  (a  typical  riometer  frequency)  as  a  fimction  of  corrected  geomagnetic  (CGM)  latitude 
(^),  corrected  geomagnetic  lon^tude  (6),  corected  geomagnetic  time  (T),  sunspot  number  (R), 
and  m<mth  (M).  This  is  subsequently  converted  to  a  value  Ai^  which  is  a  function  of  magnetic 
index  Botii  the  CCIR  (1980)  and  Hatfield  (1987)  describe  how  this  value  is,  in  turn, 
convert^  to  an  (q)o*ating  frequency  at  oblique  incidence  on  the  aurwal  D-region,  • 

Foppiano’s  model  of  mcmthly  median  abstxption  combines  the  effects  of  awOTal  substorms 
occurring  over  an  entire  m<mth  at  a  given  hour  of  the  day,  for  a  given  month  of  the  day,  and  far  a 
given  level  of  solar  activity  (specified  by  the  12-month  running-average  sunspot  number,  R). 
Although  the  median  absorption  is  implicitly  related  to  magnetic  activity,  through  its  dq>endence 
on  sunspot  number,  die  mo^l  does  not  have  an  e3q>licit  dependence  on  magnetic  activity; 
consequently,  the  large  inoeases  in  auroral  activity  tend  to  be  obscured  in  the  averaging  process. 
C<mversely,  a  few  laige  substorms  may  contaminate  the  quiet-day  predictions.  This  model  is  not 
af^qpriate  for  use  with  an  HF  ray  tracing  code  when  the  structure  of  the  refracting  regions  is 
defin^  as  a  function  of  magnetic  activity  in  a  3-hour  time  frame. 
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Accordingly,  SRI  develop^  a  model  (Vondrak  et  al..  1978)  diat  closely  follows  Fo{^iano’s 
formulation,  but  adds  an  average  magnetic-activity  dq>endence.  The  aurwal  absorptitm  on  an 
HF  ray  padi,^Aur .  is  calculated  in  three  stei»  as  follows: 

1.  Compute  Ql,  the  percentage  probability  that  one-way-vertical  absorption  at  30  MHz 
exce^  1  dB.  This  parameter  was  chosen  by  F<9piano  because  much  of  die  avail¬ 
able  auraral-absorptitm  data  is  presented  in  terms  of  Qi. 

2.  Using  an  empirical  reladonsh^)  between  Qi  and  An,  compute  An  and  ccmvmt  to  the 
related  K^i-dqimident  value,  Agp.  by  using  tl»  ii^iut  value  of  Kp. 

3.  Convert  the  (me-way-vertical  absorption  at  30  MHz  (i.e,  A^)  to  absupticm  at  die 
desired  operating  frequency  for  the  angle  of  inddmce  of  the  ray  on  the  D-region. 

Further  details  about  diis  model  can  be  obtained  from  Smith  and  Hatfield  (1987). 

To  assess  the  effects  of  spcvadic-E  on  received  signal  strength,  one  must  take  account  of  the 
fact  that  spo-adic-E  is  partially  reflecting  and  partially  transparent;  thus,  some  signal  attenuatum 
is  incurred  in  reflection  from,  o'  penetration  though,  the  layer.  The  reflection  loss  is  denoted 
hoe  as  Lj^ ,  and  the  transmissitm  (or  obscuration)  loss  by  L^.  IWo  methods  of  estimating  Lg 
and  Lg  are  provided  in  AMBCOM'  the  PhiUqw  model  (Fhillqis.  1963),  and  the  Sinno  model 
(Sinno  et  al.,  1976).  The  user  selects  the  mediod  to  be  used  by  specifying  an  input 

For  the  Phillips  method,  Lg  and  are  given  respectively  by 


Lj,  =  lOlogPa,  (dB) 

=  10  log  (1  -  Pa,)  (dB) 


(73) 


where  f^jEs  is  the  spivadic-E  blanketing  frequency,  f^Es  is  the  maximum  vertical  incidence  ordi¬ 
nary-ray  frequency  fa*  the  spmadic-E  layer,  and  the  probability  P  is  determined  by  equation  [37], 
as  is  done  in  lONCAP. 

For  the  Sinno  model,  the  empirical  formula  fcx  the  reflection  coeffidrat  r  is  given  by 


n 


1  +  lOp" 

2.077  _  _ 2.077 _ 

logp^  ”  log  fgEs  -  logf|,Es 


Pb  = 


fbBs 

fJEs 


(74) 


39 


Then  is  given  by 


Lg  =  -20logr  •  20  log  (1  ♦  lOp")  (dB)  (75) 


a  0.55  (a  typical  value  for  summer  midnight),  n  »  8.  When  p^  =0.9  (a  typical  daytime 
value),  ji  a  50,  such  that  Lg  increases  sharply  for  f  ^  >fJBa  ;  this  implies  a  ducker  layer  widi  a 
relatively  weU-definedcriti^frequaicy.  At p  ai(f^s  af^s),  s2idBforallji  .Toesti- 
mate  die  (discuratitm  loss  L„,  it  is  assumed  diat  the  transmissicm  coefficient  x  is  given  by 

t  =  ,/irrFJ  (76) 

giving  1.0  as 

Lo  *  -20  logr 
-  - 10  logll 


(dB) 


1  +  lOp' 


(77) 


Whenj]  a  8  (summer  night),  Ljt  -Lgfocp  a  0.67;  that  is,  the  layer  is  equally  reflecting  and  trans¬ 
parent  when  the  equivalent  vertical  incidence  frequent  is  about  two>thirds  of  foEs.  The  parameter 
fgEa  is  obtained  in  these  formulas  from  Leftin  et  al.  (1967).  The  values  of  f^a  are  estimated  as 
a  function  of  latitude  (lat),  day,  and  night  from  an  mq»essu»i  due  to  Kolawole  (1978) 


f^B 


0.5  +  0.2| 


\lat\  2  70 
liatl  i.  50 


65  /gEs  night 
9  fJBa  day 
UgfchSOj  ^  f^gr|70_-jiatlj  50<|iat|<70 


(78) 


Smith  and  Hatfield  (1987)  generally  i»efer  the  Sinno  model. 

The  ground  losses  are  calculated  as  in  lONCAP  at  each  reflectitm  point. 

The  field  strengdi  in  AMBCOM  is  calculated  differendy  than  in  the  othm*  analytical  pro¬ 
grams.  The  field  strength,  ,  fo'  each  mode  is  obtained  by 


Bg  -  ^l2onpo  ♦  10* 


where  the  powm*  density  is 


Po  = 


£t£t 

1 


(79) 


(80) 
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II  nH 

I  -  4n{P')^  *  *  h*  Ip*  h*  1,*  Ip.  (81) 

with 


p' 

sz 

Group  Padi  (m) 

^D2 

* 

divogence  loss  for  P2  high-^y  modes 

= 

absorpticm  loss  terms 

= 

F1-F2  deviative  absorption 

■*Aur 

ss 

auroral  absorpticm 

io 

= 

Es  obscuraticm  loss 

i* 

= 

Es  reflection  loss 

SS 

ground-reflectUm  loss 

Pc 

= 

power  transmitted  in  watts 

9t 

SS 

transmit  gain  relative  to  an  isotropic 

In  AMBCOM,  is  converted  into  decibel  units  by 

20  loge^  •  Ef  •  Pp  *  85.8  (dB>l  mv/m)  (82) 


relative  to  1-watt  radiated;  or  equivalently 

Eg  •  +  115.8  (dB>l  iiv/m) 

whoePc  is  the  power  transmitted  in  dB,  relative  to  1  watt,  and  G^  is  the  transmit  antenna  gam  in 
dB ,  relative  to  an  isotropic.  This  is  the  equation  for  field  strength  for  a  half-wave  d4>ole  over  ground 
(Nielson  et  al.,  1967).  ^  AMBCOM,  this  field  strength  is  converted  to  that  of  an  isotropic  antenna 
by  adding  8.6  dB  onto  the  free-space  loss.  In  AMBCOM,  the  field  strengths  fw  each  indi^^ual  mode 
are  not  summed  to  obtain  (me  value  rqiresenting  a  particular  frequracy  as  in  the  other  {HOgrams 
ccmsidered  here.  In  AMBCOM,  the  antenna  gaii»  are  interpolated  fi'om  tables  fweacdi  antenna  type. 
These  tabtes  are  given  at  27  take-off  angles  and  15  frequencies. 

STATISTICAL  DATA  PROCESSING 

The  CCIR  (1990b)  describes  a  {Heliminary  standardized  procedure  for  comparing  fmedicted 
and  observed  sky-wave  signal  intensities  at  ft^uencies  between  2  and  30  MHz.  Su(di  compari¬ 
sons  {H’ovide  informaticm  cm  the  accuracy  of  pr^ction  methods  and  desirable  improvements  to 
them.  To  obtain  comparable  results,  the  CCIR  suggests  diat  the  same  data  set  and  procedure  be 
used  when  making  compaiiscms  between  predicted  and  observed  field  strengths.  They  also 
suggest  that  dieir  latest  database  be  used.  Although  this  allows  the  results  to  be  compared  to 
(tiioa  by  using  their  database,  and  even  diough  most  of  the  prc^grams  compared  here  will  be  like 
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tiuwe  that  dM  COR  might  be  invested  in,  die  PENEX  database  will  be  used  instead.  They  also 
suggest  that  tbe  cmi^iarisons  should  include  the  ccamt,  the  mean  differences,  and  the  standard 
deviadmi.  The  data  screening  program  described  below  outputs  these  parameters,  as  well  as 
mhei'  parameters  useful  in  detotnining  die  accuracy  of  a  prediction  program. 

In  the  cmnpariscm  of  a  (Hx^am  against  measured  data ,  it  is  highly  desirable  to  subdivide  the 
database  into  subsets  according  to  variables  influencing  the  predicted  and  observed  results  (e.g., 
path  length,  seascm,  m<mth,  year,  gemnagnetic  latitude,  sunspot  number,  local  time  at  path 
mu^int,  etc.).  The  COR  (i990b}  suggest  diat  die  data  be  subdivuted  according  to  the  finequen- 
des  transmitted,  path  great-circle  distance,  geomagnetic  latitude,  sunspot  number,  local  time  at 
path  mu^Kunt,  caigin  of  data,  and  ratio  of  transmitted  frequency  to  (nedicted  mcmthly  median 
basic  MUF.  They  suggest  that  the  following  parameters  should  be  transferred  to  the  comparison 
program:  month,  year,  sunspot  number,  circuit  identiffer,  frequency,  great-circle  distance, 

24  predicted  hourly  values  (always  mondily  median  values)  of  slcy-wa^^  field  strength  in  dB 
relative  to  pVAn,  of  path  basic  MUF,  and  of  die  percentage  of  the  days  per  mmith  whoi  the 
frequmicy  is  below  the  path  basic  MUF.  Based  on  die  descriptions  of  the  programs  above,  the 
following  parameters  should  also  be  transfmed  to  the  comparison  |xx^am:  solar  zenith  angle, 
cosine  of  the  solar  zenith  angle.  E-layer  critical  frrequency,  B-layer  MUF,  secant  of  the  angle  of 
incidence  on  the  D-layer,  critical  frequency  of  the  ^xaradic-E  layer  fe^3 ,  the  sporadic-E  layer 
blanketing  frrequency  f^JEs ,  and  the  magnetic  index  I^,.  To  accmnplish  this,  a  computm’ 
program  called  D  ASCR3  (acronym  for  Data  Screening  3)  will  be  u^.  This  program  will  be  the 
cornerstone  of  the  PENEX  HP  {xediction  [vogram  comparistm. 

DASCR3 

DASCR3  is  a  program  designed  to  perform  data  screening  and  statistical  compariscxi  of  two 
large  matrices  of  obsovations — the  observed  data  and  the  predicted  data.  For  each  set  of 
matrices,  iqi  to  10  sets  of  information  are  read  in  on  screening  propositimis  to  be  satisfied  and 
limits  on  a  selected  variable.  In  turn,  a  portion  (ff  each  matrix  is  re^  in  and  tested  for  each  set  of 
propositicms.  For  each  subset  satisf^g  a  given  set  of  conditions,  the  variable  to  be  analyzed  is 
stoi^  tempm-arily  on  disc.  The  next  pmtion  of  each  matrix  is  then  read  in  and  screened,  and  the 
good  observaticms  are  added  to  those  already  on  disc.  When  the  entire  matrix  has  been  screened, 
data  are  then  read  into  ewe,  and  the  difference  (w  residual)  between  the  two  matrices  is  taken. 
These  arrays  are  then  swted,  to  ensure  maximum  computer  efficiency  for  the  statistical  evalua- 
tiem.  Finally,  a  statistical  evaluation  is  performed  wi  the  screened  data  and  their  residuals. 

DASCR3  has  been  used  extensively  in  past  propagation  prediction  model  analysis.  The 
following  samples  are  taken  fixnn  this  past  wwk  to  show  DASCR3  analytical  outputs  (Sailws, 
Moision  &  Brown,  1981).  An  example  of  the  output  from  DASCR3  is  shown  in  figure  4.  In  this 
sample,  the  FTSA-l  MUF  prediction  is  compared  to  observed  data.  The  preposition  to  be 
satined  is  the  data  to  be  evaluated  fw  the  month  1  (January).  In  the  printout,  the  observed  data 
are  represented  by  column  A  and  the  predicted  values  represented  by  column  B.  The  residual 
(the  observed  data  minus  the  predicted  value)  is  given  by  column  D.  The  relative  residual  is 
given  by  column  D/A,  and  the  absolute  relative  residual,  by  column  ABS  (D)/A.  The  left-half 
side  of  the  page  shows  the  statistics  calculated  for  each  of  these  columns  and  includes  the 
parameters  suggested  by  the  CCIR.  hi  addition,  the  cwrelation  coefficient  between  the  observed 
and  predicted  data  is  given.  Included  also  are  the  sl(pe,  intwcq>t,  and  mean  square  errw  of 
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linear  r^reasion.  In  this  example,  288  data  points  were  selected  frinn  4668  data  points  by  DASCR3 . 
The  average  absolute  relative  residual  for  this  case  is  25.9%.  As  this  is  an  early  applicati(Hi  of  the 
data  screening  fM-pgram,  this  particular  example  does  not  include  the  cq>ability  to  determine  a  proba¬ 
bility  distributkm  {presenting  the  residuals  described  below. 

SCREENING  DATA  BASE 

In  the  following  mcamples,  each  computer  program  was  run  to  fX'oduce  a  database  come- 
sptMKiing  to  the  observed  database.  Auxiliary  information  ou^utted  by  the  predictimi  program 
for  screening  in  this  particular  pplication  included  the  following:  universal  time  of  prpagation; 
mtmth;  year;  sunspot  number;  padi  length,  in  Idlometers;  geographic  latitude  and  longitude  of 
the  path  midpoint;  the  local  time  at  path  midpoint;  the  path  (xientation  with  respect  to  North;  the 
^lomagnetic  latitude  at  each  of  the  control  points;  the  predicted  MUF;  E-Layer  MUF;  F-layo- 
MUF;  FOT,  HPF,  and  path  identification. 

Before  the  actual  data  screening  started,  data  points  in  both  observed  and  predicted  bases 
cmresponding  to  observed  values  at  the  extremes  of  the  particular  measurement  site  were 
removed.  In  other  worxls,  obviously  pathological  outlying  data  were  removed. 

ANALYSIS  OF  RESIDUALS 

Introduction 

An  indication  of  the  accuracy  of  the  numerical  predictions  of  a  parameter,  the  MUF  in  this 
case,  can  be  obtained  from  a  study  of  the  residuals  between  observed  data  and  predicted  values. 
The  terms,  residual,  relative  residual,  and  absolute  relative  residual  are  used  with  the  following 
standard  meaning: 


Residual  =  (Observed  Datum)  -  (Predicted  Value) 

Relative  Residual  =  Resit^al - 

Observed  Datum  (g4) 

Absolute  Relative  Residual  -  Residual 

Observed  Da  turn 

Certain  statistical  measures  of  these  terms  have  proven  useful  in  past  ionospheric  studies  comparing 
predicted  and  obso^ed  data.  These  include  the  following: 

1.  The  average  residual  (avg.  res.) 

2.  Root-mean-square  residual  (rms  res.) 

3.  Mean  absolute  error  of  the  residual  (mae  res.) 

4.  Average  relative  residual  (avg.  rel.  res.) 

5.  Root-mean-square  relative  residual  (rms  rel  res.) 

6.  Mean  absolute  error  of  relative  residual  (mae  rel.  res.) 

7.  Average  absolute  relative  residual  (avg.  abs.  rel.  res.) 
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8.  Correlation  coefficient  between  obso^ed  and  predicted  values 

9.  Standard  OTW  the  estimate  of  linear  regression 

10.  The  fit  of  the  residual  distribution  to  an  empirical  probability  distribution. 

Examples  oi  each  of  these  parameters,  as  produced  by  DASCR3  (excq>t  the  last),  are  shown  in  fig¬ 
ure  4  (Sailors,  Moisitm  &  Brown,  1981). 

Description  of  Parameters 

The  average  residual  and  the  average  relative  residual  locate  the  center  of  the  distributions  of 
error,  and  they  are  sometimes  referred  to  as  die  "bias”  in  the  estimate.  Figures  5  and  6  (Sailors, 
Moision,  &  Brown,  1981)  illustrate  the  average  residual  and  the  avo-age  relative  residual, 
respectively,  as  a  fimction  of  month  fcM*  die  four  programs  compared.  In  this  example,  MINI- 
MUF-3.5  is  shown  to  have  the  smallest  bias;  who'eas,  HPMUI^4  tends  to  always  predict  high 
by  as  much  as  3.5  MHz,  (v  17.5%. 

The  mean  absolute  errors  of  the  residual  and  the  relative  residual  are  a  measure  of  the  range 
of  errm.  The  errors  are  the  first  moments  about  the  avm'age  residual  and  average  relative 
residual,  respectively.  They  fx-ovide  information  about  the  range  of  variation.  Figures  7  and  8 
(Sailos,  Moision,  &  Brown,  1981)  are  examples  of  these  two  parameters,  respectively,  for 
MINIMUF-3.5.  They  are  displayed  as  bars  about  the  average  residual  (bias)  as  a  function  of 
month;  however,  figure  5  (S^ms,  Moision,  &  Brown,  1981)  shows  that  during  the  equinox 
months  of  March  and  Sqitember,  the  range  of  variation  in  the  error  is  greater  than  the  other 
months. 

The  average  absolute  relative  residual  is  a  measure  of  the  average  magnitude  of  the  eraor. 
Figure  9  shows  a  plot  of  the  average  absolute  relative  residual  as  a  function  of  month  for  the  four 
programs  being  compared. 

The  root-mean-square  residual  and  relative  residual  are  measures  of  the  dispersion  in  the 
errm*.  In  fact,  die  RMS  residual  and  the  RMS  relative  residual  are  the  standard  deviations  of  the 
error  about  the  origin  (zero  bias),  and  diey  are  related  to  the  standard  deviation  about  the  mean 
acoxdingto 


02  =  Vj  -  V?  (85) 

where  ^2  is  the  mean  square  error  (the  square  of  the  RMS  erTw),  and^  i  is  the  bias.  When  the  bias 
is  small,  w  nearly  zero,then  the  standard  deviation  and  the  RMS  oror  are  nearly  the  same.  Other¬ 
wise,  the  rms  eirw  is  larger  than  die  standard  deviation.  Figures  10  and  11  (Sailors,  Moision,  & 
Brown,  198 1 )  are  examples  of  the  RMS  residual  and  RMS  relative  residual,  respectively,  f<x  the  four 
programs  being  compared  to  as  a  fimction  of  month.  MINIMUF-3.5  has  the  lowest  RMS  errOT, 
reaching  its  highest  value  of  4  MHz  (plus  1 2%)  during  October,  whereas,  HFMUFES4  has  its  lowest 
values  during  the  summer  months,  but  has  the  highest  RMS  error  during  the  winter  months. 

A  measure  of  the  degree  of  association,  or  the  closeness  of  fit,  between  variables  is  given  by 
the  correlation  coefficient.  It  indicates  the  strengdi  of  the  tendency  for  high  (or  low)  values  of 
one  variable  to  be  associated  with  high  (or  low)  values  of  the  other  variable.  Figure  12  (Sailcxs, 
Moision,  &  Brown,  1981)  is  an  example  of  the  corelation  coefficients  as  a  function  of  month 
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Figure  5.  Average  residual  (bias)  as  a  function  of  month. 
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Figure  6.  Avo-age  relative  residual  (relative  bias)  as  a  function  of  month. 
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Figure  7.  Average  residual  (bias)  for  MINIMUF-3.5  with  the  mean 
absolute  error  about  the  relative  residual. 
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Figure  8.  Average  relative  residual  (relative  bias)  for  MINIMUF-3.5  with 
the  mean  absolute  error  about  the  average  relative  residual. 
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Figure  9.  Magnitude  of  the  error  (average  absolute  relative  residual)  as  a  function  of  month. 
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Figure  10.  Rms  error  (in  MHz)  as  a  function  of  month. 
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Figure  11.  Rms  relative  eiT<M'  (in  percent)  as  a  function  of  month. 
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Figure  12.  Correlation  coefficients  as  a  function  of  month. 
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for  the  four  programs  being  compared.  In  this  example.  HFMUI^4  generally  has  the  highest  (xx-- 
relatitm  coefficient  with  MINIMUF-3.5  also  showing  consistently  high  values.  A  low  value  of  the 
correlation  of  an  auxiliary  parameter  may  mean  that  the  auxiliary  parameter  is  not  important  in  the 
determination  of  accuracy.  In  this  example  from  a  previous  accuracy  study,  it  was  found  that  the 
E-MUF  was  not  important  in  the  determination  of  the  MUF;  hence,  in  subsequent  accuracy  studies 
of  the  MUF.  the  E-MUF  was  not  used  as  an  auxiliary  parameter. 

A  description  of  the  nature  of  the  relationship  between  variables  is  called  regression  analysis. 
Regression  analysis  is  concerned  with  the  problem  of  describing,  or  estimating  the  value  of.  one 
variable,  called  the  dependent  variable,  on  the  basis  of  one  or  mwe  other  variables,  called 
independent  variables.  In  other  cases,  regression  may  be  used  merely  to  describe  the  relationship 
between  known  values  of  two.  or  more,  variables. 

Regression  analysis  that  involves  the  determination  of  a  linear  relationship  between  two 
variables  is  referred  to  as  simple  linear  regression.  Here,  the  variable  y  is  given  as  y  =  a  +  bx. 
where  x  is  the  independent  variable,  and  y  is  the  dependent  variable.  The  coefficients,  a  and  b. 
are  determined  in  the  regression  analysis.  A  measure  of  the  success  of  linear  regression  analysis 
is  the  standard  error  of  the  estimate  give  by 


•'y.x 


(86) 


where  is  the  standard  deviations  in  the  observed  datum,  and  the  Y  is  the  correlation  coefficient 
between  the  observed  data  and  the  predicted  values.  If  the  relationship  is  truly  linear,  then  the  bias 
of  the  estimate  should  be  removed  (or  made  nearly  zero).  An  estimate  of  the  standard  error  of  the 
mean  is 


^y.x 


(87) 


A  measure  of  the  error  in  the  regression  coefficient  is  given  by 


(88) 


where  is  the  standard  deviation  in  the  predicted  values.  The  parameter  n  is  the  p<^ulation  size. 
Figures  13  and  14  (Sailors.  Moision.  &  Brown.  1981)  show  the  standard  error  of  the  estimate  of  lin¬ 
ear  regression,  and  of  standard  errOT  of  mean  in  linear  regression,  respectively,  as  a  function  of 
month.  When  figure  13  is  compared  to  figure  12  (Sailors.  Moision.  &  Brown.  1981).  the  largest 
change  occurs  for  HFMUFES4.  Very  little  change  is  shown  for  ITSA-l  with  local  time  bq)e.  MINI- 
MUF-3.S  shows  some  changes  frxr  some  mondis.  but  not  all.  Figure  13  shows  that  linear  regression 
has  removed  much  of  the  bias  in  the  predicted  MUFs. 
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Figure  13.  Standard  error  of  estimate  of  linear  regression  as  a  funcdrm  of  month. 
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Figure  14.  Standard  errm  of  the  mean  of  linear  regression  as  a  function  of  month. 
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In  the  ^ppUcatitm  used  to  describe  the  data  screening  program,  it  turned  out  that  HFMUFES4 
had  the  laig^  bias  (7.2%),  rms  errcnr  (8.3%)  and  magnitude  of  error  (26.09%).  An  explanaticm 
was  given  for  HFMUFES4’s  poor  resets.  The  fact  that  HFMUFES4  showed  a  large  r^uctitm  in 
stamlard  errm  of  regressicm  as  a  functicm  of  nmge,  as  compared  to  die  rms  error,  led  to  the 
cooclusiiMi  that  the  model  was  in  mor;  that  is,  the  is  the  linear  portion  of  the  MUF 

model.  Had  the  results  not  shown  a  reduction  in  the  standard  error  of  regression,  as  compared  to 
rms  error,  then  die  mror  in  the  MUF  model  would  have  been  due  to  the  non-linear  portion — the 
M-factor  model.  Possibly  as  a  result  of  this  study,  newer  predictitm  programs,  such  as  lONCAP, 
did  not  use  the  set  oi  foP2  coefficients  that  were  used  in  HFMUFES4. 

Exampte  of  DASCR3  to  Improve  a  Modd 

Figures  15  through  20  describe  the  accuracy  of  QLOF  Version  2.0  as  a  function  of  local  time 
at  the  path  mic^int  by  using  several  of  the  statistical  measures  of  errm*  discussed  above  (Sailms 
and  Moisitm,  1987).  The  bias  is  nearly  zero,  from  6  to  about  17  local  mean  time  (LMT).  Since 
QLOF  is  basically  a  daytime  model,  these  results  are  encouraging;  however,  at  night,  &e  modd 
is  mwe  than  1.0  MHz  (20%)  low  on  the  average,  and  the  relative  rms  errm  is  more  than  30%,  as 
compared  to  20%  for  daytime.  At  night,  the  average  absolute  relative  residual  is  as  high  as  30%, 
as  compared  to  values  around  15%  during  the  day.  At  night,  there  are  some  hours  when  tiie 
correlation  coeffidmit  has  a  significant  decrease.  The  error  at  night  was  probably  because  the 
LUF  in  QLOF  Version  2.0  is  set  to  2  MHz  at  night. 

Because  of  this  error  at  night,  QLOF  was  modified  to  indude  nightly  D-r^on  absoqpticm  by 
using  a  modd  due  to  Wakai  (1961;  1971)  and  Wakai  et  al.  (1971).  This  was  accomplished  by 
adding  Wakai’s  nighttime  sunspot  number  dependence  into  the  existing  absoptitm  modd.  llie 
nighttime  sunspot  number  dqwndenoe  is  givra  by 

Jn»0. 025 ♦(1+0. 013 *12x2)  .  (89) 

The  new  sunspot  number  dependence  replaces  the  existing  dependeace  whenever  the  solar  zenith 
angle  at  the  control  point  exceeds  approximatdy  103  degrees.  The  QLOF  Vsrsion  2.1  showed  cmi- 
siderable  improvement  over  Vsrsion  2.0.  Table  9  ctmtains  a  comparison  of  the  accuracy  of  Vsrsitm 
2.0  and  Versitm  2.1  ova:  die  entire  database  of  observed  LOFs.  It  is  divided  according  to  tite  time 
of  day.  In  each  cell  in  die  table,  two  values  are  givoi  far  the  particular  statistical  parameter;  die  first 
value  is  for  Version  2.0;  the  second,  is  Version  2.1.  This  database  ctmsisted  of  1814  LOF  monthly 
median  LOF  obsovations  over  a  range  of  path  lengths,  seasons,  geographical  locatitms,  times  and 
sunspot  numbers.  It  is  described  in  detail  by  Sailors  and  Moision  (1987). 

Table  9.  Comparison  of  the  accuracy  of  QLOF  Vsrsions  2.0  and  2. 1  fm  day  and  night. 


Statistical  Parameter 

HmeofDay  | 

Daytime 

Nighttime 

Bias  (MHz) 

.682/.545 

1.16/.063 

RMS  Error  (MHz) 

2.02/1.94 

1.80/1.26 

Correlation  Coeffident 

.845/.854 

.716/.762 

Numbo'  of  Data  Points 

1,371  (75.6%) 

443  (24.4%) 
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Rgure  IS.  QLOP  Vsnion  2.0  average  residual  as  a  function  of  midpatfa  local  time. 
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Figure  16.  QLOF  Version  2.0  root-mean-square  residual  as  a  fimction  of  midpath  local  time. 
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Figure  17.  QLOF  Versioo  2.0  average  relative  residual  as  a  function  of 
midpadi  local  time. 


1  2  3  4  5  6  7  89Sttt213  14  1516l7BQ2021222324 


MIOPATH  LOCAL  TIME 


Figure  18.  QLOF  Version  2.0  root-mean-square  relative  residual  as  a 
fimction  of  mic^ath  local  time. 
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Figure  19.  QLOF  Vmion  2.0  average  absohite  relative  residual  as  a 
fiincticm  of  mk^ath  local  time. 
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Figure  20.  QLOF  Versitm  2.0  (xrrelaticm  confident  as  a  fimctioa  of 
mu^^atfa  local  time. 
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The  feaiUtopreaentsd  in  table  9  show  dutt  the  Version  2.1  is  markedly  reduced  from 
diat  of  Venum  2.0.  This  it  eqjecdally  noticeable  in  the  ni^ttime  results.  Similarly,  dioe  is  an 
iaqxovemeitt  in  the  rms  error  for  >^aaion  2.  The  Version  2.1  bias,  as  funcdmi  rtf  mic^Mth  local 
is  shown  in  figme  21.  This  figure  should  be  compared  to  figure  22,  which  shows  similar 
results  for  Version  2.0.  Immediately  notioeidele  is  the  reduction  in  the  bias  in  Vaai<»  2.1  for 
mir^padi  local  times  after  1700  hours  and  bdbre  0700  hoivs.  Figures  21  and  22  also  imlicate  the 
stiuidanl  deviation  of  tiie  residual  distributimi  at  each  hour. 

Enqeirical  Error  ProbabOity  Distribution 

An  additional  c^eability  in  tiie  data  scremiing  program,  which  may  (xove  useful  in  the 
FENEX  project  is  its  ability  to  fit  an  empirical  distribution  fiinctimi  to  tiie  probability  distribu¬ 
tion  r^resenting  the  residuals  (tiie  errors)  between  the  observed  parameters  and  the  correqxmd- 
iqg  predicted  parameters.  The  residuals  for  a  particular  model  are  fit  to  a  Johnsrm  syston  <tf 
freqpiency  cinrves  (Jolmsai.  1949)  by  using  an  algorithm,  which  is  credited  to  Hill  et  al.  (1976), 
Hill,  Hm.  ft  Holder  (1976),  Hill  &  Wheeler  (1981),  Dodgaon  ft  Hill  (1983),  that  uses  the 
nmtiiods  of  moments  to  (tiitain  the  required  patameters.  This  distrilmtion  rqiresents  all  univariate 
distribution  systems.  Its  sinqilicity  of  calculiuum,  once  tiie  Jrtimson  parameters  have  been 
detomined,  makes  it  adiqitable  to  miniconq;)uter  and  microcomputer  applications  because  the 
tranaformation  of  tiie  Johnami  variables  to  tiie  nonnal  system  allows  ^  use  of  normal  probabil¬ 
ity  algoritiuns  in  its  triplication.  The  nature  of  distribution  and  an  example  of  its  qiplication 
f^ow.  The  discussion  of  tiie  determinatirm  of  its  parameters  can  be  found  in  Sailors  (1987a; 
1987b). 

The  Jttimson  curves  are  an  enqiirical  family  of  curves  satisfying  the  following  choami 
conditions:  (1)  they  should  be  easy  to  evaluate  once  tiieir  parameters  are  determined;  (2)  tiiey  are 
a  moiwtonic  function  of  y,  where  y  s  x(-|)/^  Xia  a  scale  factor,  |  is  a  locatimi  factor  of  the 
dutributioi,  and  X  is  the  variable  being  rqiresented  by  tiie  distribution;  (3)  the  range  of  values 
of  f(y).  conesponding  to  tiie  actual  range  of  values  of  y,  should  be  from  -  <»  to  -koo  ,  and  (4)  the 
reaultii^  system  of  distributions  of  y  (and  so  of  x)  should  indixle  distributions  of  most,  if  not 
all,  (tf  tiie  kinds  mmountered  in  collected  data.  The  Jttfmson  syston  of  freqi^i^y  curves  ermsist 
(tf  tile  following  tiuee  types: 


the  lognormal  system  (or  5^) :  z  =  y  +  6  5<x, 

the  unbounded  system  (or  ;  z  *  y  +  6  -«»sxs« 


(90) 


the  bounded  system  (or  :  z  =  y  +  6  In 


JSiL 


[  i*X-x 


,  5<x<5+X 


where  z  is  the  standardized  normal  variate  in  eadi  case.  The  parameters  Y  and  &  determine  the 
shqie  of  the  distributkm  of  x . 
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Figure  21.  QLOF  Vsnioa  2.1  avenge  residual  and  standard  deviation 
of  the  residuals  as  a  function  of  local  time. 
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Figure  22.  QLOF  Version  2.0  avo-age  residual  and  standard  deviation 
of  the  residuals  as  a  function  of  midpath  local  time. 
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To  decide  which  one  of  the  three  Johnson  families  should  be  used  for  a  given  set  of  data,  the 
usual  procedure  is  to  obtain  the  data  estimates  of  the  third  and  fourth  statistical  mmnents  about 
the  mean — the  skewness  yPi  and  kurtosis  Hiese  are  then  plotted  on  a  graph  such  as  fig¬ 
ure  23.  Also  shown  are  other  ccmimon  sampling  distributions:  normal  (N).  Student’s  t.  Rayleigh 
(R),  and  ganuna  (w  x^)-  Figure  24  shows  data  plotted  in  the  (Pi,  P2)  plane  by  Sailca^  (1981)  that 
show  the  {X'opagatim  {H’operties  of  the  lowest  observable  frequency  (LOF)  measured  cm  the 
France-to-Iceland  path  during  1975.  In  the  case  of  the  LOF  d^  the  distribution  was  shown  to 
be  the  Sb  type  of  the  Jtdinscm  curve. 

The  data  screening  program  automatically  determines  which  one  of  the  three  Jc^mson 
families  should  be  used  by  a  given  set  of  data.  It  then  detennines  the  parameters  fm*  that 
particular  curve.  An  algcxithm  called  JNSN,  known  by  the  Royal  Statistical  Society  as  algoithm 
AS  99.  is  used  (HiU  et  al..  1976;  Hill  &  Wheeler,  1981;  Dodgson  &  Hill.  1983:  Griffiths  &  Hill. 
1985).  This  algorithm  uses  the  first  four  sample  mcnnents  to  determine  the  type  of  Johnscm  curve 
aad  its  parameters.  When  the  moments  are  luge,  the  method  of  moments  are  not  always 
statisti<^y  efficient;  consequently,  alternate  methods  are  used  to  obtain  a  second  set  of  parame¬ 
ters.  These  methods  include  (1)  maximum  likelihood  (Hahn  &  Shapiro,  1967;  Johnson,  1949), 
(2)  the  use  of  quantiles  or  pocentiles  (Hahn  &  Shapiro,  1967;  Johnson,  1949),  and  (3)  the  use  of 
frequency  moments  (Ord,  1972).  The  particular  method  used  depends  tm  the  Johnson  curve 
being  evaluated.  A  chi-square  test  of  fit  is  used  to  choose  between  the  sets  of  parameters  so 
chosen  (Hahn  &  Shapiro,  1967;  Williams,  1950).  A  5%  level  of  significance  is  used. 

After  determining  the  Johnson  distribution  parameters  by  using  the  data  screening  program, 
it  might  be  applied  in  one  of  two  ways.  The  fir^  way  is,  that  for  a  given  probability,  it  might  be 
desired  to  know  the  error  in  the  mo(^.  The  sectmd  is,  that  given  a  certain  errm*  in  the  model, 
what  is  the  corresponding  {M'obability?  The  algmidims  necessary  for  these  two  q>plications  are 
contained  in  Sailtxs  (1987a)  and  disused  below. 

In  the  first  ^)plication,  the  given  probability  is  converted  to  the  corresponding  normal 
standard  deviate  by  using  the  algorithm  functicni  PFND  (Beasley  &  Springer,  1977;  Griffith  & 
Hill,  1985).  Then,  the  corresponding  Johnson  deviates  are  found  by  using  the  algorithm  AJV 
(Dodgson  &  HiU,  1983;  Griffiths  &  Hill.  1985;  HiU,  1976;  HiU  &  Wheeler,  1981).  The  parame- 
to^  necessary  as  iiq)ut  are  ouq}utted  by  the  data  screening  program.  The  Johnson  deviates  are  the 
arar  for  the  model  being  employed. 

In  the  second  iq>pUcation,  the  given  erra-  is  converted  to  normal  standard  deviates  by  using 
the  second  algorithm  SNV  due  to  HUl  (1976),  Dodgson  &  HiU  (1983),  Griffiths  &  HUl  (1985), 
HiU  &  Wheeler  (1981).  Then  the  corresponding  probabiUty  level  can  be  found  by  using  a  ncffmal 
integral  algorithm  (Hill,  1973)  caUed  ALNORM.  This  particular  algorithm  has  the  capabiUty  to 
calculate  mther  the  uppo*  or  lower  taU  area  of  the  standardized  normal  curve  correspmiding  to 
any  given  argument 

A  sample  of  the  Johnson  distribution  appUcation  is  given  in  figure  25.  This  figure  of  the 
MUF  model  in  HFBC84,  shows  the  residual  variation  as  a  function  of  path  range.  The  predicted 
residual  is  given  fw  seven  different  standard  ncxmal  deviates  (snv)  and  their  cmresponding 
{xrdbabiUty  levels.  The  residuals  range  from  values  that  might  occur  from  0.1%  to  ^.9%  of  the 
time. 

These  same  tools  wiU  be  used  in  the  analysis  and  presentation  of  the  PENEX  field  strength 
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Figure  23.  Region  in  (^i.  ^2)  plane  fw  the  Johnson 
system  of  curves. 


Figure  24.  LOF  propagation  properties  fra*  the  France  to 
Iceland  path,  October  1975. 
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Figure  25.  Predicted  residual  for  the  HFBC84  MUF  for  the  given  standard  normal 
deviates  and  their  (xuresponding  probability  levels  with  path  range. 


CONCLUSIONS 

This  rqxMt  has  reviewed  how  seven  HF  [HXK}agation  prediction  prc^grams  derive  field 
strength.  This  review  included  three  empiiictd  programs  (Medusa  PROPHET,  FTZ,  and  FT2i4), 
and  four  analytical  programs  (HFTDA,  lONCAP,  ASAPS  and  AMBCOM).  Although  Medusa 
PROPHET  uses  the  FTZ  i^roach,  the  derivation  of  fm  and  fi  is  quite  different. 

All  of  these  analytical  iq)proaches  contain  similar  structures.  The  Erst  term  is  a  constant  that 
has  been  derived  from  summing  a  number  of  constants  along  with  several  conversion  factors.  A 
secmid  tom  accounts  for  the  frequency  that  is  in  use,  and  the  tiiird  tom  is  a  loss  tom.  There  are 
two  basic  analytical  iq^roaches  in  pr(^agation  prediction,  the  Lucas  iq)proach,  which  spawned 
versions  from  ESSA-FI^Al  throu^  lONCAP,  and  the  ray  tracing  q>proach  used  in  At^COM. 

The  prediction  schemes  produce  median  predicticms  of  the  rms  field  strength.  When  FTZ 
produces  a  value  (e.g.,  65  dB  above  1  mkrovolt  per  meter),  it  is  a  median  value  fen:  a  given  hour 
in  a  given  month.  The  median  value  is  bounded  by  tiie  description  of  the  distribution  of  the  errm 
about  that  median.  The  primary  value  in  the  PEN^  comparisons  is  in  detomining  the  size  of 
tte  standard  deviations  and  variances,  such  that  the  median  value  is  bounded. 
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The  implementation  of  the  data  screening  program  DASCR3  will  allow  the  development  and 
generation  of  a  powerful  statistical  description  of  the  measured  field  strength  characteristics  and 
of  how  accurately  the  seven  candidate  programs  predict  observations.  It  offers  all  the  statistical 
requirements  suggested  by  the  CCIR  (i990b)  for  the  determination  of  the  accuracy  of  a  field 
strength  prediction  program.  Its  ability  to  allow  and  store  up  to  40  different  auxiliary  variables 
allows  die  comparison  to  be  sub-divided  into  the  required  categories.  In  fact,  17  auxiliary 
variables  for  this  effort  were  identified  as  desirable  parameters  to  be  added  to  the  data  screening 
database.  These  variables  include:  the  month;  year;  sunspot  number;  circuit  identifter  (name  and 
path  transmitter  and  receiver  coordinates);  frequency;  great-circle  distance;  24  predicted  hourly 
values  (always  monthly  median  values)  of  sky-wave  field  strength  in  dB  relative  to  1  pV/m,  of 
path  basic  MUF,  of  the  percentage  of  the  days  per  month  when  the  fi’equency  is  below  the  path 
basic  MUF,  of  solar  zenith  angle  and  cosine  of  the  solar  zenith  angle  at  path  midpoint,  of  E-layer 
MUF;  24  predicted  hourly  values  at  each  reflection  point  (control  point)  of  E-layer  critical 
frequency,  secant  of  the  angle  of  incidence  on  the  D-layer,  critical  frequency  of  the  sporadic-E 
layer  foEs,  and  the  sporadic-E  layer  blanketing  frequency  fbEs;  and  the  eight  3-hour  magnetic 
index  values.  Useful  statistical  parameters  produced  by  DASCR3  that  can  describe  the 
accuracy  of  the  predicted  field  strength  include  the  following;  average  residual  (bias);  root- 
mean-square  residual  (standard  deviation);  average  relative  residual  (relative  bias);  root-mean- 
square  relative  residual;  average  absolute  relative  residual  (magnitude  of  the  error  in  the  model); 
correlation  coefficient  between  observed  and  predicted  values;  standard  error  of  the  estimate  of 
linear  regression;  the  constants  necessary  to  represent  the  residual  distribution  by  a  Johnson 
probability  distribution,  and  its  corresponding  test  of  fit  information.  DASCR3  usage  will  also 
allow  the  determination  of  possible  improvements  that  might  be  made  to  these  field  strength 
prediction  programs. 
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